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Abstract 

The lace expansion is a powerful perturbative technique to analyze the critical behavior 
of random spatial processes such as the self-avoiding walk, percolation and lattice trees and 
animals. The non-backtracking lace expansion (NoBLE) is a modification that allows us to 
improve its applicability in the nearest-neighbor setting on the Z'^-lattice for percolation, 
lattice trees and lattice animals. 

The NoBLE gives rise to a recursive formula that we study in this paper at a general 
level. We state assumptions that guarantee that the solution of this recursive formula 
satisfies the infrared bound. In two related papers, we show that these conditions are 
satisfied for percolation in d > 11, for lattice trees in d > 16 and for lattice animals in 
d > 18. 
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1 Introduction and results 

1.1 Motivation 

We use the non-hacktracking lace expansion (NoBLE) to prove the infrared bound for several 
spatial models. The infrared bound implies mean-field behavior for these models. The classical 
lace expansion is a perturbative technique that can be used to show that the two-point function 
of a model is a perturbation of the two-point function of simple random walk (SRW). This 
result was used to prove mean-field behavior for self-avoiding walk (SAW) [261128] . percolation 
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[Ml [29] , lattice trees and lattice animals [25] , oriented percolation [Mj , the contact process 
[3I1E3I and the Ising model [42] in high dimensions. 

Being a perturbative method in nature, applications of the lace expansion typically ne¬ 
cessitate a small parameter. This small parameter is often the inverse of the degree of the 
underlying graph. There are two possible approaches to obtain a small parameter. The first 
is to work in a so-called spread-out model, where long-, but hnite-range, connections over a 
distance L are possible, and we take L large. This approach has the advantage that the results 
hold, for L sufficiently large, all the way down to the critical dimension of the corresponding 
model. The second approach applies to the simplest and most often studied nearest-neighbor 
version of the model. For the nearest-neighbor model, the degree of a vertex is 2d, which has 
to be taken large in order to prove mean-field results. 

For the self-avoiding walk (SAW) on the nearest-neighbor lattice, Kara and Slade (1991) 
[261IM] proved the seminal result that dimension 5 is small enough for the lace expansion to 
be applied, thus that mean-field behavior holds in dimension d > 5. This result is optimal 
in the sense that we do not expect mean-field behavior of SAW in dimension 4 and smaller. 
The dimension 4 thus acts as the upper critical dimension. Results in this direction, proving 
explicit logarithmic corrections, can be found in a series of papers by Brydges and Slade (some 
also with Bauerschmidt), see [11] and the references therein. 

For percolation, we expect mean-field behavior for dimensions d > 6. Kara and Slade also 
proved this result down to d > 6 for the spread-out model with sufficiently large L [241129] . For 
the nearest-neighbor setting, Kara and Slade computed that dimension 19 is large enough. 
These computations were never published. Through private communication with Takashi 
Hara, the authors learned that in a recent rework of the analysis and implementation the 
result was further improved to d > 15 for percolation. 

To obtain the mean-field result also in smaller dimensions above the upper critical di¬ 
mensions, we rely on the NoBLE. In the NoBLE, we explicitly take the interaction due to 
the last edge used into account. Doing this, we reduce the size of the involved perturbation 
drastically and are able to show the mean-field behavior for dimensions closer to the upper 
critical dimension. 

In this paper, we formalize a number of assumptions on the general model and prove that 
under these assumptions the two-point function obeys the infrared bound. The derivation 
of the model-dependent NoBLE and the verification of the assumptions are not part of this 
article. We use the generalized analysis to obtain mean-field behavior for the following models: 
lattice trees in d > 16 and lattice animals in d > 18 HZ], and percolation in d > 11 [T8] . 

A NoBLE analysis consists of four steps: Eirstly, for a given model, we derive the pertur¬ 
bative lace expansion. Secondly, we prove diagrammatic bounds on the perturbation. Thirdly, 
we analyze the expansion to conclude the infrared bound given certain assumptions on the 
expansion. In our analysis, we derive diagrammatic bounds on the lace-expansion coefficients 
in terms of simple random walk integrals, which can be computed explicitly. This allows us 
to compute numerical bounds on the lace-expansion coefficients. The fourth step consists of 
the numerical computation of these SRW-integrals. 

In the accompanying papers HZl [18], we perform the hrst two steps for percolation, as 
well as for lattice trees and lattice animals. In this paper and in a model-independent way, we 
perform the analysis in the third step and explain the numerical computations of the fourth 
step. The numerical computations and the explicit checks of the sufficient conditions for the 
NoBLE analysis to be successful are done in Mathematica notebooks that are available on 
the website of Robert Fitzner [14] . 
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The analysis presented in this paper is an enhancement of the analyses performed by 
Hara and Slade in m (see also [Ml ESI ESI ES] for related work by Hara and Slade), and 
by Heydenreich, the second author and Sakai in [32] • This paper is organized as follows: 
In Section 11.21 we first introduce simple random walk (SRW) and non-backtracking walk 
(NEW). In Section fl.4l we state the two basic NoBLE relations that are perturbed versions 
of relations describing the NEW. Then, we state the results for percolation, lattice trees and 
lattice animals proved in [laiiB]. 

In Section [2l we first explain the idea of the proof at a heuristic level. Then, we state all 
assumptions required to perform the analysis in the generalized setting and state the result we 
prove in this document, namely the infrared bound. We close in Section [2.51 with a discussion 
of our approach. 

In Section [3l we prove the technical cornerstone of the analysis, namely, that we can 
perform the so-called bootstrap argument. For the analysis in Sections 1^1151 we use a simplified 
NoBLE form of the two-point function that allows us to present the analysis in a clearer way. 
Thus, we also state the assumptions in Section [2.21 in terms of this simplified characterization. 

In Section S[ we explain how to derive the simplified NoBLE equation starting from the 
NoBLE equation for the two-point function, and reformulate the assumptions of Section 12.21 
into assumptions on the NoBLE coefficients that are derived and bounded in the accompa¬ 
nying papers [nuBj. 

Section [5] is devoted to the numerical part of the computer-assisted proof. We explain 
how we compute bounds on the required SRW-integrals. In Section [5.31 we explain the ideas 
to bound the NoBLE coefficients that are used for all models that we consider. We end this 
paper with a general discussion. 

1.2 Random walks 

We begin by introducing the random walk models that we perturb around and fix our notation. 

1.2.1 Simple random walk 

Simple random walk (SRW) is one of the simplest stochastic processes imaginable and has 
proven to be useful in countless applications. For a review of SRW and related models, we 


refer the reader to [SSlEZlIlll. 


An n-step nearest-neighbor simple random walk on is an ordered (n -|- l)-tuple oj = 
{iOQ,u]i,uj 2 , ■ ■ ■ ,ujn), with u!i € Z'^ and \\uji — cjj+iHi = 1, where ||x||i = Unless 

stated otherwise, we take ujq = 0 = (0,0,..., 0). The step distribution of SRW is given by 



( 1 . 1 ) 


where 6 is the Kronecker delta. For two functions f,g: Z'^ i—M and n € N, we define the 
convolution f -k g and the n-fold convolution /*"■ by 



( 1 . 2 ) 


and 


/-(x) = (/*(-!) * /)(x) = (/*/*/*...* /)(x). 


(1.3) 
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We define Pn{x) as the number of n-step SRWs with ujn = x, so that, for n > 1, 

Pn{x) = ^ 2dD{y)pn-i{x -y) = 2d{D -kpn-i){x) = (2d)”R)*”(x). (1.4) 


We analyze this function using Fourier theory. For an absolutely summable function /, we 
define the Fourier transform of / by 

fik) = f{x)e^^"' for k € [-vr,7r]'^, (1.5) 

x£Z‘^ 


where k ■ x = kiXi, with inverse 

fix) = 


' {—7r,7r)^ 


( 1 . 6 ) 


We use the letter k exclusively to denote values in the Fourier dual space (—tTjTt)'^. We note 
that the Fourier transform of /*"'(x) is given by fik)^ and conclude that 


1 

Pn(fe) = (2d)"'i)"'(A:), with F>(fc) = - 'Y, cos(fcJ. 

t=i 

The SRW two-point function is given by the generating function of pn, i.e., for z G C, 

OO ^ 

Czix) = ^p„(x) 2 ;”, and C^{k) = - - 

n=0 ^ ~ ‘^d.ZlJ[k) 

in x-space and A:-space, respectively. We denote the SRW susceptibility by 

x"^'^(^) = 4(0) = ^ 


(1.7) 


( 1 . 8 ) 


I-2dz 


(1.9) 


with critical point Zc = l/(2(i). By the form of Cz{k) in (11.81) and using that 1 —cos(t) ^ t^/2 
for small t € R, we see that Cz^ik) = [1 — D{k)]~^ 2(i/||/c||2 for small k, where || ■ ||2 denotes 

the Euclidean norm. Since small k correspond to large wave lengths, the above asymptotics 
is sometimes called the infrared asymptotics. The main aim in this paper is to formulate 
general conditions under which the infrared bound is valid for general spatial models. 


1.2.2 Non-backtracking walk 

If an n-step SRW ui satisfies Ui a;j _|_2 for all z = 0,1,2,... , n — 2, then we call w non¬ 
hacktracking. In order to analyze non-backtracking walk (NBW), we derive an equation 
similar to dUD. The same equation does not hold for NBW as it neglects the condition that 
the walk does not revisit the origin after the second step. For this reason, we introduce the 
condition that a walk should not go in a certain direction b in its first step. 

We exclusively use the Greek letters i and n for values in {—d, —d -|- 1,..., —1,1, 2,... , d} 
and denote the unit vector in direction z by G Z'^, e.g. {ef)i = sign(z)d|^| j. 
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Let bn{x) be the number of n-step NBWs with loq = 0,ujn = x. Further, let b‘^{x) be the 
number of n-step NBWs lo with Un = x and coi 7 ^ e^. Summing over the direction of the first 
step we obtain, for n > 1 , 


bnix)= 5^_i(x-be,). (1.10) 

Further, we distinguish between the case that the walk visits —in the first step or not to 
obtain, for n > 1 , 


bn{x) = b^‘'{x) + b‘^_i{x + e,). ( 1 . 11 ) 

The NB W two-point functions and B^^ are defined as the generating functions of bn and 
bn, respectively, i.e.. 


CXD 00 

Bz{x) = Y bn{x)z^, B^(x) = Y (1.12) 

n=0 n=0 

Using (II.IOD and (|1.11|) for the two-point functions gives 

Bz{x)=5o,x + z Y K{x + e,), B^{x) = B^^x)-\-zB'^^{xe,). (1.13) 

Taking the Fourier transform, we obtain 

B,{k)=l+ Y 13^{k) = B-\k) + ze-'^’^^B‘-^{k). (1.14) 


In this paper, we use C^'^-valued and x C^'^-valued functions. For a clear distinction 
between scalar-, vector- and matrix-valued quantities, we always write valued functions 
with a vector arrow (e.g. v) and matrix-valued functions with bold capital letters (e.g. M). 
We do not use {1,2,..., 2d} as index set for the elements of a vector or a matrix, but use 
{—d, —d-|- 1 ,..., — 1 , 1 , 2 ,... , d} instead. Further, for a A: € (— 7 r, 7 r)'^ and negative index 
i € {—d, —d -|- 1 ,..., —1}, we write k^ = —k\i\- 

We denote the identity matrix by I € and the all-one vector by 1 = ( 1 , 1, ... , 1)^ € 

Moreover, we define the matrices J,D(A:) € by 

(J)t,K = dt_K and (D(fc)),,K = ( 1 . 15 ) 

We define the vector Bz{k) with entries {Bz{k))i, = B}(k) and rewrite ( 11 . 141 ) as 

B^{k) =1 + zl^t){-k)Bz{k), Bz{k)l = 3Bz{k) + zt){-k)Bz{k). ( 1 . 16 ) 

We use JJ = I and D(A:)D(—/c) = I to modify the second equation as follows: 


Bz{k)l = 3t){k)t){-k)Bzik) + z33t){-k)Bzik) = 3{'D{k) + z3)t){-k)Bzik) (1.17) 


which implies that 


t)i-k)Bzik) = Bzik) i){k) + z3 


-1 


Jl. 


(1.18) 
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We use J1 = 1 and then combine (jl.lSI) with the first equation in (ll.lGp to obtain 

1 


B,{k) = 


l-zP 


1 -1 


1 


Then, we use that 


D(A:) + zJ 


1 -1 


1 — z 


D{k) + z3 

t){-k)-z3] , 


and r^D(—A:)l = 2dD{k) to conclude that 


B,{k) = 


1 


1-^2 


l_2dz^^^ '\- +{2d-l)z‘^-2dzD{k) 

1 — 2^ 


(1.19) 


( 1 . 20 ) 


( 1 . 21 ) 


The NEW susceptibility is x'^®'^(z) = Bz{Q) with critical point Zc = l/(2d — 1). The NEW 
and SRW two-point functions are related by 

1-^2 1 


B,{k) = 


1 + (2d - 1).2 1 _ 


1-^2 


C_ 


1 -I- (2d — 1)^2 i+(2d-i)z' 


Xk), 


so that 




( 1 . 22 ) 


(1.23) 


2d-ll-Zl(A:) 

This link allows us to compute values for the NEW two-point function in x- and A:-space. A 
detailed analysis of the NEW, based on such ideas, can be found in [T6] . 


1.3 General setting 

We consider general models defined on the d-dimensional hypercubic lattice TX. For these 
models, the two-point function i—>■ M is the central quantity. The two-point function 

is defined for parameters z G [0, ^^c) where Zc acts as the critical value. As for the SRW and 
NEW, the susceptibility (52(0) diverges as z approaches Zc from below. The behavior of Gz 
and Gz as z Zc is of special interest. We use the NoELE to prove that the two-point 
function of the general model is a small perturbation of the critical NEW two-point function 
(|1.23p and thereby obeys the infrared bound. 

To do this, we define as the two-point function of the model where is being avoided. 
The precise definition depends on the model. For NEW, is avoided in the first step, for 

percolation is the two-point function of the model defined on the graph Z'^ \ {e^}. For the 
NEW, these two-point functions are linked by the relations (ll.lip . In the NoELE, we adapt 
these two relations for the general model with model-dependent perturbations, and bound 
the arising perturbation coeffcients. For d > 2, the NoELE gives rise to functions 
and for l,k^ ■ ■ ■ j all mapping from Z'^ to M, and a function Hz ■ K+ —>■ IR+, 
such that, for all x G Z*^ and 2 G [0, Zc), 

Gz{x) = So^x + ^zix)fiz{^o,y ^z{y))G‘'^{x — ye,^), (1-24) 

Gz{x) = G‘'^{x) + fizG^Xx — e ^)^ n(,’''(y)G^(x — y-|-e^)-I-^'"(x). (1-25) 
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In our applications, the variable /Xj, is closely related to the main parameter z, but it is not 
equal. Therefore, in the analysis we use a second parameter p,z that allows us to control the 
critical value in the models under consideration. For example, for percolation we use 

fLp = p, fJ-p = P^piei not connected to 0 | the bond (0,ei) is vacant). (1-26) 

(See Section fl.4.11 where the percolation model is formally introduced.) 

Our goal is to understand the behavior of Gz, where we consider the functions pz, 
as given. Applying the Fourier transformation on ()1.24l) and ()1.25l) gives 

Gz{k) = l + Ezik) + pz Yl il + K{k)>~'"^G^zik), (1.27) 

te{±i,...,±d} 

Gz{k) = G‘,{k)+pze^^^G-^{k)+ Y f^‘'z^{k)e-'^’"^G^,{k)+ E^{k). (1.28) 


We define the vectors Gz{k),E{k) and T(A:) and the matrix tiz{k) by 

{Gz{k))^ = G‘',{k). {^{k))^ = K{k), {E{k))^ = E^,{k). {tlz{k))^^^ = tlT{k)- ( 1 . 29 ) 

Then, we can rewrite ()1.28p using vectors and matrices as 

Gz{h)l = ^z{k) + pzt>{k)3^z{h)+tlz{k)t^{-k)^z{k) + ^{k). (1.30) 


We obtain from this that 


1 -1 


Gz{h)=J^{k) D(A:)+/XzJ + n,(A;) {Gz{k)l - E{h)). 


(1.31) 


That this matrix inverse is well defined will be shown in Section Ol Next, we rewrite (11.2711 
in vector-matrix notation and solve for Gz{k) as 


Gz{k) =1 + %{k) + pz{l + ^z{k)ft>{-k)Gz{k) 


with 


-1 


=l + Ez{k) + pz{l + ^z{k)Y T>{k) + pzi + Ilz{k) {Gz{k)l-Ez{k)) 


1 -1 


^ ^z{k) — PzY-zik)) T){k) + pz^ kkl.z{k) ^zik) 


i-Pz{i + i>z{k)Y D(fc) + M + n,(fc) 


-1 


1 


^z{k) 

l-Fz{k)' 


1 -1 


^z{k) 1 + ^z{k) — PzY- + ^z{k)) D(^) + A^zJ + nz(^) ^z{k), 


Fz{k) ■.=pz{l + ^z{k))'^ t){k) + pzJ + Uz{k) 


-1 


1 . 


(1.32) 


(1.33) 

(1.34) 


When comparing (jl.32ll - (ll.34ll to its equivalent for NEW in (jl.lOp . we see that (jl.32ll reduces 

to (|1.19l) when taking pz = z, Ez{k) = (T^(/c))k = (nz(/c))t,K = {Ez{k))i^ = 0 for all z,k,i,K. 
Our analysis is based on the intuition that the NoBLE coefficients are small in high dimensions. 
The majority of our work is to quantify this statement. 




































Rewrite of the two-point function. We use another characterization of the two-point 
function Gz{k) to perform the analysis. We extract all contributions involving constants and 
D{k), by defining Cp^z, c^,z,ct<s>,z, otp^z, R^,zik), RF,zik) such that 

^zik) := c^^z + a,i.^zD{k) + R^^zik), (1.35) 

Fz{k) := Cp^z + aF,zD{k) + Rp^z{k). (1.36) 


Then, 


Gz{k) 


__ 

1 - Fz{0) + Fz{0) - Fz{k) 

T Oi^^zF{k) F R<i,z{k) 

^*2(0)/Gz(0) + otp^z[^ — D{k)\ + Rp^z{0) — RF,z{k) 


(1.37) 


In Section m we show how we transform (ll.33p - (ll.34p into (ll.35l) - (ll.36p . Up to that point, 
we only work with the representation (ll.35l) - (|1.36l) as it simplifies and shortens the analysis. 
The quantity in (| 1.3711 reduces to the NBW-equivalent (11.211) . when we set a^^z = RF,z{k) = 
Rt>,z{k) = 0 and c^^z = 1 — z^-,OiF,z = “^dz, ^> 2 ( 0 )/G^( 0 ) = 1 -|- {2d — l)z^ — 2dz. 


1.4 Results for specific models 

In this section, we describe the results that our method allows to prove. These results are 
proved in two accompanying papers mi EH]. 


1.4.1 Percolation 

Percolation is a central model in statistical physics and is a very active field of research since its 
rigorous definition by Broadbent and Hammersley in 1957 m, where this model was proposed 
to describe the spread of a fluid through a medium. General references for percolation are 
[6l 1201136] . A review of recent results can be found in [211131j and the references therein. We 
consider Bernoulli percolation on the hypercubic lattice. We use the definition of [431 Section 
9]: To each nearest-neighbor bond {x,y} we associate an independent Bernoulli random 
variable which takes the value 1 with probability p and the value 0 with probability 

1 — p, where p G [0,1]. If n^x,y} = 1; then we say that the bond {x, y} is open, and otherwise 
we say that it is closed. A configuration is a realization of the random variables of all bonds. 
The joint probability distribution is denoted by Pp with corresponding expectation Ep. 

We say that x and y are connected, denoted by x i —>■ y, when there exists a path 
consisting of open bonds connecting x and y, or when x = y. We denote by ^{x) the random 
set of vertices connected to x and denote its cardinality by |^(x)|. The two-point function 
Tp{x) is the probability that 0 and x are connected, i.e., 

Tp{x) = Pp(0 —> x). (1.38) 

By translation invariance Pp(x ^—> y) = Tp{x — y) for all x, y € We define the percolation 
susceptibility, or expected cluster size, by 

Xip) = = ^P • ( 1 - 39 ) 
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We say that the system percolates when there exists a cluster ^(x) such that |^(x)| = oo. 
We define 6{p) as the probability that the origin is part of an infinite cluster, i.e., 

e{p) = ¥p{\^iO)\ = oo). (1.40) 

For d> 2, there exists a critical value Pc € (0,1) such that 

Pc{d) = inf{p I e{p) > 0}. (1.41) 

Menshikov in 1986 [39], as well as Aizenmann and Barsky in 1987 [2|, have proven that the 
critical value can alternatively be characterized as 

Pc{d) = sup {p I x{p) < oo} • (1-42) 

The percolation probability p i—)• 9{p) is clearly continuous on [t),pc), and it is also continuous 
(and even infinitely differentiable) on (pc, 1] by the results of [3] (for infinite differentiability 
of p !->■ 9{p) for p € (Pc 5 1], see [lQ|). Thus, the continuity of p i-)- 9{p) on is equivalent to 
the statement that 9{pc{d)) = 0. 


Critical exponents. We introduce three critical exponents for percolation. It is widely 
believed that the following limits exist in all dimensions: 


7 = — lim 
/3 = — lim 
1/6 = 


logx(p) 

p/'Pc log(lp-Pcl)’ 

log0{p) 

p\/c log(lp-Pcl)’ 

logPp,(|'^( 0 )| > n) 


— lim 

n—>oo 


logn 


(1.43) 

(1.44) 

(1.45) 


A strong form of (|1.43ll . (11.441) and (|1.45p is that there exist constants c^, cg, C 5 G (0, 00 ) such 
that 


X(p) = (l + o(l))c^(pc-p) ^ asp/'pc, (1.46) 

Gip) = {I + o{l))co{p - pcf asp\pc, (1.47) 

Ppc(l^(0)l > ^) = (1 + o{l))csn~^^^ as n ^ 00 , (1-48) 

and is expected to hold in all dimensions, expect for the critical dimension d = dc, where 


logarithmic corrections are predicted. The constants c-^,ce and cs depend on the dimension. 
We say that these exponents exist in the bounded-ratio sense when the asymptotics is replaced 
with upper and lower bounds with different positive constants. Further, it is believed that 
there exist p and ci, C 2 such that 

Tp.(x) = (1 + 0 ( 1 ))^^^, fp,(A:) = (l + o(l))y^^, (1.49) 

where ci and C 2 depend on the dimension only. For percolation, the existence of many more 
exponents is conjectured and partially also proven. See m Section 2.2] for more details. Our 
main result for percolation is formulated in the following theorem: 
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Theorem 1.1 (Infrared bound for percolation). The infrared bound fp^{k) < A 2 {d)/[1 — D{k)] 
for some constant A 2 {d) holds for nearest-neighbor percolation in dimension d satisfying d > 
11 . As a result, the critical exponents 'y,/3,5 and rj exist in the bounded-ratio sense and take 
their mean-field values 7 = /3 = l,(5 = 2 and rj = 0. 

Theorem II.II is proved by combining the model-independent results proved in this paper, 
with the mo del-dependent results as proved in [IS]. There, we also state and prove related 
results on percolation, such as the existence of the so-called incipient infinite cluster and the 
existence of one-arm critical exponents. Further, we derive numerical upper bounds on the 
critical percolation probability Pc{d). 

The critical exponents for percolation have received considerable attention in the litera¬ 
ture. For the critical exponents, it is known that /3 < 1 and 7 > 1 for all d > 2, see Chayes 
and Chayes m and Aizenman and Newman [3]. Further, we know that if fd and 7 exists, 
then fd € (0,1] and 7 G [1, 00 ), see [20l Section 10.2 and 10.4]. In high dimensions, we expect 
mean-field behavior for percolation. Namely, we expect that for all dimensions d > 6 , the 
critical exponents correspond to the exponents of the regular tree given by 7 = /3 = 1 , and 
7 = 0. (For the definition of 7 for percolation on trees, see Grimmett [201 Section 10.1].) 
Alternatively, we can interpret 7 = /3 = 1, and 7 = 0 as the critical exponents for branching 
random walk, see the discussion in [31] . An important step to prove mean-held behavior 
for percolation is the result of Aizenman and Newman [3] that the hniteness of the triangle 
diagram, dehned by 


^(Pc) = (rp,*rp^*TpJ(0), (1.50) 

implies that 7 = 1 . This triangle condition also implies that (3 = 1, see [3]. In particular, this 
implies that p 9{p) is continuous. 

Kara and Slade [23] use the lace expansion to prove that 7 = 0 in Fourier space as well as 
the hniteness of triangle diagram for d > 7 in the spread-out setting with a sufficiently large 
parameter L. In the spread-out setting all bonds {x,y} with \x — y\ < L are independently 
open or closed. This is an optimal result in the sense that mean-held behavior is not expected 
in d < 6 , see [36], where Toulouse argues that the upper critical dimensions dc, above which 
we can expect mean-held behavior, equals 6 . For mathematical arguments why dc = 6 , see 
Chayes and Chayes m, Tasaki |45j or [311 Section 11.3.3]. 

For the nearest-neighbor setting, Kara and Slade proved mean-held behavior in sufficiently 
high dimensions [23|. Later, they numerically verihed that d = 19 is sufficiently high by 
adapting the proof of the seminal result that self-avoiding walk in dimensions d > 5 satishes 
the infrared bound. In private communication with Takashi Kara, the authors have learned 
that in a recent improvement of their numerical methods, the mean-held result was established 
for d > 15, and thus this implies Theorem ll.il The proof for both these results (d > 15 and 
d > 19) were never published. 

Let us briehy explain how percolation hts into our general framework. For percolation, in 
[18] . we perform the non-bracktracking lace expansion (NoBLE) for the two-point function 
Tp{x) and Tp(x) = Pp(0 ^—> x without using ef). Further, we bound the coefficients arising in 
this expansion and check that all general assumptions used in the present paper are satished. 
In the NoBLE, we further identify that p,p = pPp(e € "^(0) | {0, e} is vacant). For our analysis 
we also require a bound on p.p = p (recall ()1.26l) ). 
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1.4.2 Lattice trees and animals 

A nearest-neighbor lattice tree (LT) on is a finite, connected set of nearest-neighbor bonds 
containing no cycles (closed loops). A nearest-neighbor lattiee animal (LA) on is a finite, 
connected set of nearest-neighbor bonds, which may or may not contain cycles. Although a 
tree/animal A is defined as a set of bonds, we write x € A, for x G Z'^, to denote that x is an 
element of a bond of A. The number of bonds in A is denoted by |A|. We define tn\x) and 
tn\x) to be the number of LAs and LTs, respectively, that consist of exactly n bonds and 
contains the origin and x G Z'^. We study LA and LT using the one-point funetion gz and 
the two-point function Gz defined as 


gP = GPiO) = 

A: ASO 

gP = GP(O) = 

E 

T: T30 

(1.51) 

00 

E J"' 

00 

<5f’(U = E‘.‘.‘’W^” = 

E 

(1.52) 


n=0 A: A^0,x n=0 T: T90,x 


where we sum over lattice animals A and trees T, respectively. For technical reasons, we 
perform the analysis for the normalised two-point function Gz{x) = Gz{x)l9z- This is not 
necessary, but simplifies our analysis in the general framework and improves the numerical 
performance of our method. 

We define the LA and LT susceptibilies by 

X(“>(^) = X^^\z) = Gf\0), (1.53) 

and denote the radii of convergence of these sums by and Zc \ respectively. As for SRW 
and NEW, l/zc describes the exponential growth of the number of LTs/LAs as the number 
of bonds n grows. When we drop the superscript (a) or (t), we speak about LTs and LAs 
simultaneously. The typical length scale of a lattice tree/animal of size n is characterized by 
the average radius of gyration given by 

Rn = — - ||x|| 2 tn(a^)- (1-54) 

Critical exponents. The asymptotic behavior of tn and Gz can be described using critical 
exponents. We define three of these critical exponents for LA and LT. In doing so we drop 
the superscripts (a) and (f) as the following holds for LA and LT. It is believed that there 
exist 7 , (5, z/, ?7 and Ai, A 2 , A 3 , A 4 > 0 such that 

Xiz) = (1 + o{l))-- —, Rn = {l + o{l))A 2 -n^, (1.55) 

(1 - z/zcP 

GzPx) = + and G^^(/c) = (1-b o(l))A 4 ||/c||^"^ (1.56) 

\\x\\2 

as z zA Zc, n —>■ 00 , ||x ||2 —>■ 00 and k ^ 0, respectively. The exponents are believed to be 
universal, in the sense that they do not depend on the detailed lattice structure (as long as 
the lattice is non-degenerate and symmetric). In particular, it is believed that the values of 
7 , 5, rj and n are the same in the nearest-neighbor setting that we consider here, and in the 
spread-out setting. The constants Aj do depend on the lattice structure. 
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For LT/LA, it is believed that the critical exponents take their mean-held values above 
their upper critical dimension dc, which are 7 = 1/2, u = 1/4, rj = 0. These values correspond 
to the mean-held model of LT/LA, studied in [7]. It is conjectured in [38] that the upper 
critical dimension of LT and LA is dc = 8 . This conjecture is supported by [SO], where it 
is shown that if the “square diagram” is hnite at the critical point, as is believed for d > 8, 
then the critical exponent 7 satishes 7 < 1/2. In [9], it has been proven that 7 > 1/2 in all 
dimensions. 

For the nearest-neighbor setting that we consider, Kara and Slade give a rigorous proof of 
mean-held behavior for LT and LA in sufficiently high dimensions, see [25|. What sufficiently 
high dimensions means was not made precise. The authors learned through private commu¬ 
nication with Takashi Kara that it was not investigated in which dimension the classical lace 
expansion starts to works. In particular, Hara and Slade expected it to only be successful in 
dimensions much larger than dc = 8 . In the spread-out setting with L large enough, Hara 
and Slade in |25| proved the mean-held behavior for LT and LA in all dimension d > 8. Our 
main result for LTs and LAs is formulated in the following theorem: 

Theorem 1.2 (Infrared bound for LTs and LAs). The infrared bound Gz^{k) < A{d)/\\. — 
D{k)\ holds for some A{d) for nearest-neighbor lattice trees in dimension d satisfying d > 16, 
and for nearest-neighbor lattice animals in dimension d satisfying d > 18. As a result, 7 takes 
its mean-field value 7 = 1/2. The critical exponent t] exists in the bounded-ratio sense and 
takes its mean-field value p = 0. 

From our analysis we show directly that r/ = 0 in Fourier space, see the right side of (|1.56|1 . 
and can easily deduce that 7 = 1/2. The proof that that 7 = 0 in x-space is non-trivial and 
is explained in m using the results by Hara [22| . 

For LTs and LAs, in we perform the non-bracktracking lace expansion (NoBLE) on 
the two-point function Gz{x) and G''fix), which is the two-point function in which the LT and 
LA partially avoid e^,. In the expansion, we identify yiz = 25 ^ and fiz = zQz, where Qz = Gz(0) 
and = ^ 2 ( 0 ) 3-^6 the one-point functions. 

The expansion proves two relations for Gz{x) and G‘'^{x) that are perturbations of (I1.14p . 
The NoBLE can further be used to obtain bounds on the NoBLE coefficients and to verifies 
that the assumptions formulated in this paper indeed hold in the dimensions in Theorem ll.21 

2 Main results 

The main result of this paper is the infrared bound in a general setting. In this section, we 
first explain the idea of the proof. Then, we state the assumptions on the general model and 
prove the infrared bound under these assumptions. We close this section with a discussion of 
our results. 

2.1 Overview 

The NoBLE writes Gz{k) as a perturbation of the NBW two-point function, see ()1.32D . where 
the perturbation is described by certain NoBLE-coefficients that we denote by and 

Hz^. In the accompanying papers, we derive the NoBLE and its coefficients and prove that 
they can be bounded by a combination of simple diagrams. When we can bound these simple 
diagrams, then we are able to bound the perturbation, and thereby also to derive asymptotics 
for the two-point function. 
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Thus, we would like to bound simple diagrams for all 2 : < Zc- It turn out there exists 
a zi £ (0, Zc) such that we can bound the simple diagrams in terms of NBW diagrams 
for all z < zj. For example, for percolation, zj = l/{2d — 1). To obtain bounds also for 
2 € {zj,Zc), we use a bootstrap argument, which is a common tool in lace-expansion proofs, 
see e.g. [laiaiES]. We next explain how such a bootstrap argument works. 

We use the following minor modification of the classical bootstrap argument: 

Lemma 2.1 (Bootstrap argument). For i = 1,2,3, let z 1 — fi{z) be continuous functions 
on the interval \zi,Zc)- Further, let 7 j,rj G W be such that 1 < 7 i < Tj and fi(zi) < 7 *. If 
for z € {zi,Zc) the condition fi{z) < Tj for all i G {1,2,3} implies that fi{z) < 7 * for all 
i € (1, 2, 3}, then in fact fi{z) < 7 * for all z € [zj, Zc) and i € {1, 2,3}. 

Proof. We consider the continuous function 2 ; i-)- maxj=i^ 2,3 fi{z)/Ti and see that the lemma 
follows directly from the intermediate value theorem for continuous functions. □ 

To apply the bootstrap argument, we use three different functions that we need in order 
to bound the lace-expansion coefficients: 

(a) fi to bound the critical values of fiz, fiz and 2 ; G {zi,Zc)‘, 

(b) /2 to bound the two-point function Gz{x) in Fourier space; 

(c) /s to bound so-called weighted diagrams, such as weighted bubbles or triangles. 

We define and explain these functions at the end of this section. We first continue our 
explanation at a more heuristic level. 

For z = zj, we bound the simple diagrams using the NBW diagrams and use these bounds 
to prove that fi{zi) < 7 *. This ‘initializes’ the bootstrap argument. For 2 : G {zj,Zc), we use 
the idea depicted in Figured] . Firstly, we assume that fi{z) < Fj and use this assumption 
to conclude bounds on various diagrams consisting of combinations of two-point functions. 
Then, we use these bounds to bound the lace-expansion coefficients. This, in turn, enables 
us to conclude bounds on the bootstrap functions fi{z). //these bounds turn out to be 
such that fi{z) < < Fj, then we can use Lemma 12.11 to conclude that fi{z) < 7 ^ for all 
z < Zc- These bounds then imply the infrared bound. We extend the result to Zc using a 
left-continuity property of the two-point function and the NoBLE coefficients, which we prove 
in the accompanying papers as these arguments are model-dependent. 

The structure of the proof, shown in Figure dl is the reason that we are not able to 
prove the infrared bound for all dimensions above the upper critical dimensions dc- The 
perturbation, identified by the NoBLE, is for dimensions close to dc quite large and reduces 
as we increase the dimension. Thus, for high dimensions such as d > 100, it is relatively 
simple to show that the statement that fi{z) < Fj implies that /j( 2 ;) < 7 j < Fj. However, it is 
very difficult to prove such a statement for dimension closer to the upper critical dimension dc- 
The dimensions stated in Theorem ll.ll and ll.2l do not show anything specific about the model, 
but only the limitation of our technique. In fact, for example for percolation, d > dc = Q is 
the proper condition. 

Using exhaustive bounds on the model-dependent NoBLE coefficients and a more tedious 
computer-assisted proof might allow to prove the infrared bound in dimensions above, yet 
closer to, dc- Thereby, it is not clear whether it can be used to obtain the infrared bound for 
all dimension above dc for percolation, LT and LA. 

We next explain the idea of our proof in more detail, so as to further highlight the ideas 
in this paper. We continue by defining and discussing the bootstrap functions. 
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Assume a bound 



Figure 1: Structure of the bootstrap argument: Improvement of bounds 
Bootstrap functions. For the bootstrap, we use the following functions: 


fi{z) 

:= max{(2d-l)^ 2 ,c^( 2 d-l)/i^}, 

(2.1) 

f2{z) 

\Gz{k)\ 2d-l 

:= sup = sup [1 D{k)] \Gz[k)\, 

k£{ — 7V,7r)^ ^ — 

(2.2) 

hiz) 

suPxGsEy \\y\\lGz{y){GT * D*'^){x - y) 

\= max - - - , 

{n,l,s}es Cn^l,S 

(2.3) 


where > 1 and Cn^i^s > 0 ^ire some well-chosen constants and S is some finite set of indices. 
Let us now start to discuss the choice of these functions. 

The functions fi and /s can been seen as the combinations of multiple functions. We 
group these functions together as they play a similar role and are analyzed in the same way. 
We do not expect that the values of the bounds on the individual functions constituting /i 
and /s are comparable. This is the reason that we introduce the constants and Cn,i^s- 
The value of n is model-dependent. For SAW, we would use only n = 0. For percolation 
we use n = 0,1 and n = 0,1, 2 for LT and LA. This can intuitively be understood as follows. 
By the x-space asymptotics in (I1.49P and (11.561) . and the fact that (f -k f)(x) ~ ||x|| 2 ~'^ when 
d> 4 and f{x) ~ ||x|| 2 ~'^, we have that ||y|||G^(y) ~ (Gz kGz){y)- As a result, this suggests 
that 

\\ygGz{y){GT * D*^){x - y) ^ ^{Gz k Gz){y){GT * D*^){x - y) (2.4) 
y y 

so that finiteness of WygGz{y){GY * D*^){x — y) is related to finiteness of the bubble 

when n = 0, of the triangle when n = 1 and of the square when n = 2. 

The choices of point-sets 5 G 5 improve the numerical accuracy of the method. For 
example, we obtain much better estimates in the case when x = 0, since this leads to closed 
diagrams, than for x ^ 0. For x being a neighbor of the origin, we can use symmetry to 
improve our bounds significantly. To obtain the infrared bound for percolation in d > 11 we 
use 

5 = {{0,0, A}, {1,0, A}, {1,1, A}, {1,2, A}, {1,3, A}, {1,4, {0}}}, 
with A = {x G Z'^: ||x ||2 > 1}. This turns out to be sufficient for our main results. 
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2.2 Assumptions 

In this section, we state the assumptions that we need to perform the general NoBLE analysis. 
The assumptions are in terms of the simplified form of the NoBLE in (jl.37p . In Section 01 we 
derive this simplified NoBLE form of the NoBLE and translate the following assumption in 
terms of the NoBLE-coefficients. We begin with two assumptions on the two-point function 
that are completely independent of the expansion: 

Assumption 2.2 (Bound for the initial value). There exists a zj £ [0,Zc) sueh that 

2d-2 

Gz{x) < Bi/(2d_i)(x) = ^^_^ Cxi2d{x) (2.5) 

for all X £7^ and z £ [0, zj]. 

To control the growth of the two-point function as we approach the critical value Zc, we use 
the following two assumptions: 

Assumption 2.3 (Growth of the two-point function). For every x £ the two-point 
functions z Gz{x) and z !->■ Gi,(x) are non-decreasing and differentiable in z £ (0,Zc)- For 
all e > 0 there exists a constant Cg > 0 such that for all z £ (0, Zc — s) and x £ 7'^ \ {0}, 

-—Gz{x) < C£{Gz * D-k Gz){x) and therefore —G^(0) < CsGziO)'^- (2-6) 

CL/C CL/O 

For all z £ {0,Zc), there exists a constant K{z) < oo such that ll®ll 2 G^^(®) < Fi{z). 

Assumption 2.4 (Continuity). For z £ [0,Zc), z ^ fiz ond z ^ giz are continuous. 

We only consider models where the two-point function has the following set of symmetries: 

Definition 2.5 (Total rotational symmetry). We denote by Vd the set of all permutations of 
{1,2,..., d}. For V £ Vd, d £ (—1,1}'^ and x £ we define p{x; n, 6) £ Z*^ to be the vector 
with entries {p{x]n,d))j = djXyy We say that a function / : Z'^ !->■ M is totally rotationally 
symmetric when f{x) = f{p{x; n,6)) for all n £ Vd and 6 £ (—1,1}'^. 

Assumption 2.6 (Symmetry). We assume that x Gz{x),x i-)- Rp^z{x) and x !->■ Rq,^z{x) 
are totally rotationally symmetric. Further, we assume that the lace-expansion coefficients 
satisfy 


T^o) = ^'."(0), j;n^'’"(o) = j;nr'(o) (2.7) 


for all L,K £ (±1, ±2,... , id} and z < Zc. 

The following is the central assumption to perform the bootstrap. We assume that if fi{z), 
f 2 {z), fsiz) are bounded for a given 2: £ [0,Zc), then the functions ap^z,cx.s>^z, Rf,z, R<s>,z obey 
certain diagrammatic bounds. The form of these bounds is delicate and depends sensitively 
on the precise model under consideration. 

Assumption 2.7 (Diagrammatic bounds). Let ri,r 2 ,r 3 > 0. Assume that z £ {zi,Zc) is 
such that fi{z) < Tj holds for all i £ {1,2,3}. Then, Gz{k) > 0 for all k £ (—7r,7r)'^, and the 
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following bounds hold with /?, depending only on ri,r 2 ,r 3 ,(i and the model: 


(a) There exist fd,, > 

> 0, such that 




(2.8) 


e 

VI 

(2.9) 

(b) There exist > 0? such that 



^nr(x)<:dn^. 


(2.10) 


X 


(c) There exist I5 r,f, ldR,.t., fdAR.F, ^ > 0 such that 


\Rf,z{x) 1 ^Pr,f, 

X 

y ^ i.^*,2(3;)i < Pr,<s>j 

rr 

(2.11) 

'y ^ 1 'PPar,<s>j 

X 

y ] ll^lli1 Pz Par,fi (2-12) 

T 

Rf,z{0) — RF,z{k) > — ~ 

b{k)], 

(2.13) 

for all k E (—vr, tt)'^. Further, we assume that B — B >0 and B 

•' ' ’ _ ' ’ L—a.F L—AR.F ilc,« 

If AssumMion fW2l holds, then the bounds stated above also hold for z = zj, 

“ /3|q,4>| “ > 0. 

where in this case 


the constants /3, only depend on the dimension d and the model. 

To extend the infrared bound io z = Zc, we use the following assumption: 

Assumption 2.8 (Growth at the critical point). We assume that, if the bounds stated in 
Assumption fS. 7| hold uniformly for z E [zi,Zc), then z !->■ Gz{k) is left-continuous at z = Zc 
for any k 0, and that the bounds stated in Assumption \2.7\ also hold for z = Zc- 

Assumptions 12.6112.81 depend on the NoBLE and are stated in terms of its simplified form 
(|l.c{7p . In Section m we replace these assumptions by assumptions on the NoBLE-coefficients. 
We have chosen to use the form (ll.37p for the analysis, as it simplifies the presentation of the 
analysis considerably. 

2.3 Main result: Infrared bound 

To successfully apply the bootstrap argument, we require that we can improve the bound on 
the bootstrap functions. This is the content of the following condition: 

Definition 2.9 (Sufficient condition for the improvement of bounds). For 7 ,r E and 
z E [zi,Zc), we say that P{'y,r,z) holds when fi{z) < Tj for i E {1,2,3} and the following 
conditions hold: 


0 < 7i < Tj 


7 i > max 


{/i(z/),max{/3^,C/,} 


for i = 1,2,3, 

1 + /3n‘ 


1_2^/3 

^ 2d-lP. 


72 > 


2d — 1 + /3|a,$i + 


/3 -/3 

—a,F —AR,F 


2d-2 

and 73 is larger than the maximum of the right-hand sides of (|3.3ip and ()3.87p . 


(2.14) 

(2.15) 

(2.16) 
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The last condition states that the initial condition < 73 holds and that the im¬ 

provement of bounds succeeds for We do not give a formal statement at this point, as 
it is involved and would require notation that has not yet been introduced. If the bootstrap 
succeeds, then we are able to prove our main result: 


Theorem 2.10 (Infrared bound). Let k G [—vr, vr]'^, z/ G [0,Zc) and'y,T G If Assumptions 
\2.S^2.8\ and P{'y,T,z) hold for all z G [zi,Zc), then, for all z G [zi,Zc], 


with 


G,{k)[i-b{k)] < 

Gbk) < 


2d-2 
2d-I 


72 , 


m 

i/GM + [^-D{k)y 


A{d) 


min {/3 ^ 


/^c,* + /3|c.,*| + 



(2.17) 

(2.18) 

(2.19) 


We postpone the discussion of Theorem I2.1UI to Section 12.51 We continue to discuss the 
strategy of proof. 


2.4 Proof subject to a successful bootstrap 

The central statement, that we prove in Section [3l is that we can apply the bootstrap argument 
when P{'y,T,z) holds. We now formalize this statement in Proposition 12.111 

Proposition 2.11 (A successful bootstrap). Let 7 , T G M^. If Assumvtions I2.2H2.SI and 
P{'y,T,z) hold for all z G [zi,Zc), then the functions /i,/ 2,/3 defined in (I2.ip - (|2.3I) are 
continuous, fi{zj) < 7 ^ for i G {1,2,3} hold, and fi{z) < Tj for all i G {1,2,3} implies that 
fi{z) < li for all i G {1,2, 3}. 

We prove Proposition 12 . Ill in Section [3] one function /j at a time. Now we prove our main 
result. Theorem 12.101 assuming that Proposition 12.111 holds: 


Proof of Theorem \2.10\ subject to Proposition \2.11\. By Lemma 12.11 

^Gz{k)[l - b{k)] < f 2 {z) < 72 for all z G {zi,Zc). 


2d-2 


( 2 . 20 ) 


By Assumption 12.81 2 ; 1 —>■ Gz{k) is left-continuous at z = Zc for k ^ 0. From this, we conclude 
that (12.20p also holds for z = Zc and k ^ 0, which proves (|2.171) . To prove (12.18p . we use the 
bounds of Assumption 12.71 on Gz{k) as given in (11.3711 


Gz{k) = 


C$,2 T cx^^z-bi^k') -(- R^^zi^k') 




< 


< 


/^C,# P\cc,^\ “1“ I^R,^ 


^c,4> 


Gz(0) 




/^c,$ A /3|a,$| + 


( 2 . 21 ) 


1/G2(0) + [1 - b{k)] ’ 

which implies (12.1811 and derives the expression for A{d) in (|2.1911 . □ 
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2.5 Discussion 


In general, a proof using the NoBLE consists of four parts, see also Figure 12.51 (a) the 
derivation of the non-backtracking lace expansion or NoBLE; (b) diagrammatic bounds on the 
NoBLE coefficients; (c) the analysis of the NoBLE equation; and (d) a numerical verification 
of the conditions in P( 7 ,r, 2 ;), using a computer-assisted proof. 


Expansion 

General relation ^ 

Analysis/ Bootstrap 

Model dependent 

Model independent 

in accompanying papers 


in this paper 


Coefficients to 

Prior bounds 


describe the 

on the two-point 


perturbation 

function 






Bound on the 
perturbation 


Diagrammatic bonnds 


Nnmerical bounds 

Model dependent 
in accompanying papers 

Model dependent 

in accompanying Mathematica notebooks 

Bounds m 


form of diagrams 


T 


Numerical values 


Numerical computation 

Model dependent 

in accompanying Mathematica notebooks 


Figure 2: Structure of the non-backtracking lace expansion. 


Parts (a) and (b) are performed in the model-dependent papers [THl [IT]. Part (c) is 
performed here in a generalized setting. Part (d) is explained in Section (5] and is numerically 
performed in three Mathematica notebooks. In the first notebook, we compute SRW-integrals 
for a given dimension, see Sections 15.1115.21 In the second notebook, we implement bounds on 
the simplified rewrite (| 1.371) and the bound necessary for the improvement of /a, see Appendix 
El and Section [3.3.51 respectively. These two parts are completely model independent. In the 
third notebook, we use the values of the SRW-integrals and the bootstrap assumptions to 
compute numerical bounds on the diagrammatic bounds on the NoBLE coefficients. These 
bounds are then used to verify the conditions P{'^,T,z), which, when successful, imply that 
the analysis here yields the infrared bounds in the specific dimension under consideration. 
Since the bounds are monotone in the dimension, the bounds then also follow for all dimensions 
larger than that specific dimension. 

In the thesis of the first author m, the analysis was performed in two ways. The first 
was based on the x-space approach, as originally worked out by Kara and Slade in EZj. This 
approach was used by Kara and Slade in [ST] [28| to obtain that mean-field behavior holds for 
self-avoiding walk (SAW) in all dimensions d > 5. This is optimal in the sense that mean-field 
behavior for SAW in d < 4 should not be true. See m and references therein for results in 
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this direction. Further, Kara and Slade adapted their method to percolation |24j . which led 
to the famous result that mean-field behavior for percolation holds for d > 19. 

The second analysis was based on the trigonometric approach, first used for finite tori in 
[5] and worked out for in [321143j . However, it was never verified above which dimension 
this technique can be applied. Thus, it was initially not obvious to us which method would 
be numerically optimal. It was only by implementing both methods that we discovered that 
the x-space approach, combined with the NoBLE analysis, is numerically superior. This 
is the reason that we only describe this method here. In conclusion, our method is, after 
the derivation of the NoBLE, heavily inspired by that of Hara and Slade in m- We have 
benefitted tremendously from their work, as well as from the many discussion that we have had 
with Takashi Hara and Gordon Slade over the past years. From private communication, we 
have learned that Takashi Hara also managed to prove that percolation in dimension d > 15 
obeys the infrared bound, although this result has not appeared in print. We hope that our 
method, as well as the accompanying Mathematica notebooks that are publicly available [13] 
to anyone who is interested, increase the transparency of the proof of the infrared bound for 
all the models involved. 

The main difference between our work and the work by Hara and Slade in [27j is that in 
our method, the loops creating the perturbative terms are made to consist of at least 4 bonds, 
while in the classical lace expansion they could consist of immediate reversals (2 bonds). This 
makes the perturbation considerably smaller, and allows for an analysis that is to a much 
larger extent model-independent as the analysis in [27| and its adaptation to percolation. It 
also explains why our method gives reasonable results for lattice trees and lattice animals, 
models that previously had not been attempted by Hara and Slade. In discussions with 
Takashi Hara, we have found that our bounds on e.g. the triangle diagram are slightly better 
in dimension 15, whereas he has a much more sophisticated and model-dependent analysis of 
the lace-expansion coefficients. 

For the SAW, we also derived a NoBLE and implemented the bootstrap. In this way, we 
can show that mean-field behavior holds for SAW in d > 7, see m and [131 • While the proof 
for d > 7 is relatively simple, we expect that an extension of the technique to d = 5,6 will not 
produce a substantially simpler proof than that of Hara and Slade, that is already optimal in 
the sense that it proves the result in all dimension above the SAW-upper critical dimension 
4. Thus, we have not attempted to improve upon our result. 

Let us dwell a bit on the distinction between the x-space approach and the /c-space or 
trigonometric approach. We require bounds on weighted diagrams, alike 

-Rf,2(0) — Rp^zik) = — cos(A: • x)] < [1 — D{k)]p. (2.22) 

X 

We can either bound the underlying diagram directly in Eourier space or use the following 
lemma to translate it into the x-space approach: 

Lemma 2.12 (Eourier transforms and step distribntions). For a summable, non-negative 
function g that is totally rotationally symmetric, as defined in Definition \2.fA the following 
bound holds: 


g{x)[l - cos{k ■ x)] < [1 - L»(/c)]^ 5 (x)||x||^. (2.23) 
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To distribute the weight ||x ||2 over a large diagram into the weights of parts of the diagram, 
we use the relation that, for x* € 

=X] 11^*112+ (2-24) 

2=1 2=1 2=2 J = 1 2=1 2=1 

The diagrammatic bounds for the fe-space approach are very similar, due the following anal¬ 
ogous result, which is of independent interest: 

Lemma 2.13 (Split of cosines). Let t € M and tj G M for i = 1,..., J such that 
t = Then, 


.1 J i-i 

1 — cos(t) < ^[1 — cos(tj)] -|- ^ sin(fj) sin (2.25) 

2=1 2 = 2 j = l 

J 

1 — cos(t) < J ^[1 - cos(tj)]. (2.26) 

2=1 

The inequality (12.261) with a factor 2J -|- 1 is commonly used in the lace-expansion liter¬ 
ature. While reviewing the proof, the authors found that a minor adaptation improves the 
leading factor to be J . The proof of Lemmas 12.121 and 12.131 can be found in Appendix [Bj 

Let us close this section by proposing some extensions of our work. We do not manage 
to prove the infrared bound all the way down to the upper critical dimension for percolation. 
For this, we need even better arguments. One might hope that this can be done by a more 
careful analysis that compares the interacting models with memory-m walk for large values 
of m. Here, a SRW is called a memory-m walk when it has no loops of length at most m. 
Thus, NEW is memory-2. We can easily derive such a memory-m expansion for SAW, for 
percolation, LT and LA this is already more involved. However, the analysis required for this 
expansion is much more involved, and we have not tried this more general approach. One 
particular problem is that we do not know what the memory-m Green’s function is, so that it 
is harder to explicitly expand around this. We think that also for this approach a numerical 
(and thus computer-assisted) proof is necessary. 

3 Verification of the bootstrap conditions 

In this section, we prove Proposition 12.111 one function /j at a time. 

3.1 Conditions for fi 

In this section, we prove that the properties of /i in Proposition 12.111 hold. 

By Assumption 12.41 z fiz and z faz are continuous, so that z fi{^) is also 

continuous. From (j2.15j) . we conclude that fi{zj) < 71. To show that for all 2; G { zi , Zc ), 
fi{z) < Fj for all i G {1,2,3} implies that fi{z) < 71 , we prove a relation between Bf^{0) 
and Gz{0). We use the abbreviations ifz = T''(0) and where the choice of 

i G {±1,... , id} is by Assumption 12.61 not relevant. 
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Lemma 3.1 (Link between NEW and general susceptibility). Let Assumvtion \2.6\ hold and 
define 


(1 + 

1 + TT^ - ^ 2 ^’^ ’ 


SO that 


1 + '^2 

(1 + '!/’ z )/^2 + fi’z 


(3.1) 


Then, Bxfi0)^z{0) = 62 ( 0 ) for all z < Zc- 

Proof. Since D(0) = I and ^(0) = the two-point function Gz{k) in the form of (11.3211 
simplifies for /c = 0 to 


62 ( 0 ) = 


^*.( 0 ) 


1 - ^ 2(1 + ^-( 0 ))^ [d( 0 ) + ^ 2 J + n 2 ( 0 ) 

^>.( 0 ) 


1 -1 


(3.2) 


1 “ 1^2(1 + ' 02 ) 1 ^ I + + 1^2(0) 


1 -1 


Due to the simple form of I and J and the symmetry of stated in Assumption 12.61 the 
sum of each column and the sum of each row of I -|- ^z^ + 1 ^ 2 ( 0 ) equals 1 + ^z + t^^■ Thus, the 
one vector 1 is an eigenvector of I -|- ^z^ + n 2 ( 0 ) corresponding to the eigenvalue 1 + ^z + t^'' 
and we can compute that 


62 ( 0 ) 


^*.( 0 ) 


1 - ^lz{l + fiz) 


2d 

l+IMz+Tfi. 


and 


^a(O) 


2dX • 
1+X 


Solving Ea^( 0 )$ 2 ( 0 ) = G 2 ( 0 ) for A 2 for the first equality, and for Hz for the second, gives the 
desired result. □ 


The above identification allows us to improve the bound on /i as required in the bootstrap 
analysis: 

Lemma 3.2 (Improvement of /i). Let z G {zi,Zc) and 7 , T G M^. If Assumvtions \2.b^2.T\ 
and condition P( 7 , r, 2 ;) hold, then fi{z) < 71 . 

Proof. Recall that fi{z) = max{(2(i — l)fizi c^{2d — We select A 2 as in Lemma [3TT] and 

note that A 2 < {2d — 1)“^ if z < Zc as G 2 ( 0 ) = LiA 2 ( 0 )$^( 0 ) < 00 . Hence, we can compute 
that 


^ ( 2 d-l) LemmaO (1 + K){2d - 1 )A 2 M {1 + f3n.){2d - l)Xz 


1 + V’2(l + A 2 ) 


1 - + ^ 2 ) 


(2d-LA2<l _l + /3^. 


- 1 _ 

2q!— 

Using from Assumption 12.71 we obtain 


flzi‘^d-1) < fi. 


1 -|- 


1 _ 

^ 2d-lE. 


Thus, 


fi{z) <max|^^ ^ 


1 -I- fijjL 


zd Q ' ^ 1 ^d P 

which is by (12.1511 smaller than 71 when condition P{'y,T,z) holds. 


■}. 


(3.3) 


(3.4) 


(3.5) 

□ 
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3.2 Conditions for /2 


In this section, we prove that the properties of /2 in Proposition 12. Ill hold. We start with the 
required continuity. 


Lemma 3.3 (Continuity of / 2 ). The function z i->- f 2 {z) defined in (12.2p is continuous for z 
in [0,Zc). 

Proof. We follow the proof of [321 Lemma 5.3]. To show that /2 is continuous on [0,Zc), we 
prove that it is continuous on the closed interval [0,Zc — s] for any e > 0. Using Assumption 
Ea we know that for any k and z G [0, Zc — e], 


^Gz{k) 

dz 




< cfiGM? < 


(3.6) 


where we can interchange differentiation and summation as the sum is bounded in absolute 
value, as just shown. From this, we conclude that the derivative of f 2 (z) is uniformly bounded 
on [ 0 , Zc — e], which implies the continuity of /2 on [ 0 , Zc — e], □ 


We continue to prove the bootstrap for f 2 '. 

Lemma 3.4 (Improvement of / 2 ). Let z G \zi,Zc) he such that Assumvtions lg.hl|g. 7| and 
P{'j,r,z) hold. Then f 2 {z) < 72- 

Proof. Recall that f 2 iz) = sup;i,g(_j,. — D{k)]Gz{k). As already used in the proof of 
Theorem 12.101 we know that Assumption 12.71 implies that 




(3.7) 


and 


1 — Fz{k) — 4>^(0)(jr5;(0) ^ + cxp^z[f — T){k)] + Rp^z{ft) — Rp^z{k) 




(3.8) 


where we use in the last step that 1 — Fz{0) = (l> 2 ( 0 )G^ ^(0)) = > 0 by Lemma IRTI 

We conclude from this that 


\Gzik)\[l-Dik)] 


^z(fc)|[l — L){k))] ^ /^cl. + + fi\RM 

1 - 7(t) ^ 


(3.9) 


If condition P( 7 ,r, 2 ;) holds, then this is smaller than 72 , see (|2.16p . which completes the 
proof. □ 


3.3 Conditions for f^ 

In this section, we show that the function /s, defined in (12.311 . satisfies the conditions of 
the Bootstrap Lemma fLemma 12.111 . Namely, we prove that z 1 —>■ f^iz) is continuous, that 
f‘s{zi) < 73 , and that for all z G {zj, Zc), ffiz) < R for all i G {1, 2,3} implies that fsiz) < 73 . 
As this is more elaborate for f^ than for fi and / 2 , we divide the proof into multiple steps. 

The techniques of this section are an adaptation of those used by Kara and Slade to prove 
the mean-field behavior for SAW in d > 5, see [28] . The central idea needed for the adaptation 
was developed in discussions with Takashi Kara. 
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3.3.1 Rewrite of /s 

We analyze the function 


(3.10) 


nfix) = ||!,||^o,(!,)(or *o*')(i - y). 

y 


and conclude the desired results on /a using that, by definition (|2.3p . 


h{z) 


max 

{n,l,S}GS 


Cn,l,S 


(3.11) 


We bound Tiz’^ using the continuous Laplace operator A. For a differentiable function g and 
s e {1, 2,..., d}, let dsg{k) = ^g{k) and Ag{k) = Then, 

d^,Gz{k) = = - x2G,(fe)e‘'=-", (3.12) 

xeZ'i xezd- 

AGz{k) = Yl \\x\\2Gz{x)e^^'A (3.13) 

Thus, we can bound '1-Ci'\x) using the Fourier representation 

K’Kx) = /(_,,,).(- A Gz{k))&{k)G-M^-'^^-^^- (3.14) 

If we replace Gz in (I3.14p by C'i/ 2 d (recall (II.8M . then we can compute the value directly, 
see Section [3.3.31 To obtain a bound for 2: G { zi , Zc ), in Sections I3.3.4IIT331 we extract a 
dominant SRW-like contribution from Gz, which we compute directly and then we bound the 
remainder terms separately. The bounds are expressed using several SRW-integrals that can 
also be computed numerically as we explain in Section [SJ 


3.3.2 Continuity of /a 

Lemma 3.5 (Continuity). The function z !->■ fsiz) as defined in (|3.10l) is continuous in 
z G [zi,Zc). 

Proof. We fix an e > 0 and prove that {TCf’\x)) is an equicontinuous family of functions 
and is uniformly bounded for all x,n,l for all 2 ; G [0,2;c — e). This allows us to obtain the 
continuity of 2 t—> svL^^^gTCf'\x) for all sets S directly from the Arzela-Ascoli Theorem. This 
implies the continuity of 2 1 —>■ f-i{z) as the index set S over which we take the maximum in 
(ESI) is finite. By Assumption 12.31 there exists a constant K{zc — e) < 00 such that 

'Y \\x\\lGz,-eix) < K{zc-e). (3.15) 

Further, Gz^-fiO) = xi^c — e) < 00 , so that, uniformly for 2 G [ 2 /, 2c — e], 

nfi^ix) < sup(Gf * D*%)K{zc -e)< x{zc - eTK{zc - e). (3.16) 

y 
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By Assumption 12.31 


< c. E, WvWliG, G,)iy){GT * D*^){x - y) 

+nCeY.y \\y\\lGM{Gt^^KD<^+^)){x-y). (3.17) 

We use that Hw + x + y \\2 < 3(||iu||2 + ||x ||2 + Hylli) all w,x,y G Z'^ to obtain 

^ MliG, * H * G,){y) < ||u;||iG,(^n)(Zl * G,){y - w) 

yezd- y,^ 

+ 3YiG.^G,)iy-w)Diw)\\w\\l 

y,w 

+ 3YiGz * D){w)G^{y - w)\\y - w\\l 

y,w 

< 6K{zc - e)G.,-e(0) + 3G,,_e(0)2. (3.18) 

We conclude that 

< Cs{6K{zc -e) + 3G,,_,(0) + nK{zc - e))G',,_,(0)"+i < oo. (3.19) 
By the uniformity of this bound in x, we conclude that {'Hz’\x)) is equicontinuous. □ 


3.3.3 Bound for the initial point f^{zi) 

In this section, we prove that fsizi) < 73 . By Assumption 12.21 we can bound Gzj{x) < 
i^Ef^ 1 /( 2 ^)(®)- Thus, we can bound fsizj) using only SRW-quantities. We start by com¬ 
puting the derivatives of C'i/ 2 d(fc) = C{k) = [1 — D{k)]~^ and D{k), where dg denotes the 
derivative w.r.t. kg'- 


and 


We use 


a (^(u\ ^ dgDjk) ^ 2 Ef=iSia(fc^) 

s' ^[i-T)(fe)P [i-m? ’ 

.AdsD{k)?\ 


(3.20) 

(3.21) 


AD{k) 

d 

Y{dsD{k)f 

s=l 


Y^lbik) = -i^cos(fe,) = -b{k), 

S=1 

(3.22) 

d 

^2 ^sin 2 (/c,) := Z)"‘"(A:). 

(3.23) 


S = 1 


siii"(fc) = -1 (e'*- - e-“-)" = 1 - le"''- - ie'^- = 1 [1 - cos(2fc,)] (3.24) 
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to compute that 




d 



L 


- 6(2^1. 


(3.25) 


We define M{k) = D{k) — 2D^^'^{k)C{k) and conclude from the computations above that 


AC(A:) = -C{kfM{k) = -D{k)C{kf + \c{kf - ^ C{kf. 

d 2d^ 

L 

We use this representation to compute the SRW analogue of 'Hz’\x) as 

^\\y\\lCiy)iD*Ucn(^-y) 

y 

d'^k 

U\k)C"'^\k)M[k)e~^'''-^ ■ 

^ (—7r,7r)^ 


(3.26) 


&{k)C^+\k)M{k)e-^’^---^ 

(27rj“ 


' (—7r,7r)' 


&{k)C^+‘^{k) (^D{k) - ^C{k) + ^ 


d^k 

{2ttY 


(3.27) 


(3.28) 


As we explain in more detail in Section [5.11 we can numerically compute the SRW-integral 


4,Kx) := (Z1 *'*Ci7(2,))(x) = 


&{k) 


-ik-x 


d'^k 

{27rY 


'(—tTjTt)'^ [1 — £l(/c)]”' 

We use this integral to compute (I3.28p numerically as 

( |3.28P = In+ 2 ,l+l{x) - ^In+sA^) + ^'^^n+3,lA + 2eJ := Jn,l{^)- 


Thus, fsizj) is bounded by 


2d-2 snpx€sJn,iA) 

h[zi) < ——- max -, 

2tt 1 {n7,5'}E«S ^n,l,S 

By the assumption in P( 7 , T, z), this is smaller than 73 . 


(3.29) 


(3.30) 


(3.31) 


Remark 3.6 (Close to the upper critical dimension). The bound (|3.3U|1 can only be used in 
dimension d > 2(n + 3) + 1 as it uses In+ 3 ,i+i{x), which is only finite in these dimensions. 
This restricts the analysis shown here to dimensions d > dc + 3, e.g. for percolation we can 
only use this bound for d > 9 as we require a bound on 'HI’^(x) to successfully apply the 
bootstrap argument. This problem can be avoided using a different bound for the integral. 
For example, using the bound 

^ Ell - cos"(t.)] < ^ Ell - cos(fc,)| < ?[1 - D(k)] (3.32) 

S = 1 S=1 

in (I3.27p . we obtain that 

Y.WvWlCiyAD*^ * CnA -y)< In+ 2 ,l+l{x) + ^Kn+ 2 ,lA), (3.33) 

y 


where we introduce the SRW-integral Kn+ 2 ,i{x) in (13.361) below. This and other bounds 
applicable in d = dc -|- 1, dc -|- 2 , perform numerically worse than the bound in (13.311) . As 
we are not able to prove mean-field behavior in dimension dc -|- 1 , dc -|- 2 anyway, we use the 
numerically better bound (13.311) instead. 
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3.3.4 Preparations for the improvement of bounds for 

We wish to prove for all z G {zi,Zc) that fi{z) < Fj for all i € {1,2,3} implies fsiz) < 73 . 
This is the most technical part of our proof. We do this by deriving a bound on 

r fl^h 

«7(I)=/ (-AG.(t))D'(t)G;(tKi‘''^. uni 

J ( —TTjTt)*^ 

For this bound, we extract the SRW-like contributions from (— A Gz{k)) by decomposing it 
into five terms Hi ,..., We compute the SRW-like contributions Hi as in the preceding 
section and bound the remainder terms using Assumption 12 . 71 and certain SRW-integrals that 
we define next. 


SRW-integrals. Here we introduce several SRW integrals that we use to bound H'^'\x). 
Flsing the terminology introduced in Definition 12.51 we define 

E (3.34) 

’ '^G'Pd 56{-l,l}d 


Performing the sum over 5 gives cosines, so that D^^\k) is real. The following SRW-integrals 
are adaptations of the integrals used in 1271 Section 1.6]: For x € and n,/ G N, we let 



/ 

J (—77,77)' 


/ 

J (—77,77)' 

Tn,l{x) = 

/ 

J (—77,77)' 


/ 

J (—77,77)' 


b{kfC{kYb^^\k) 


d^k 


\b{k)\^c{kY\b^^\k)\ 


d'^k 




\b\k)\C{kY\b^'^\k)\\M{k)\ 


d^k 

{2ttY' 

d^k 


id’ 


(3.35) 

(3.36) 

(3.37) 

(3.38) 


where C{k) = C'i/ 2 d(fc) is the critical SRW two-point function and M{k) is defined above 
(13.261) . For any function / such that f{x) = f{p{x-,v,6)) for all v,5 (see Definition 12.51) . we 
see that 


' ( — 77 , 77 )*^ 


/(fc)e 


—\k-x 


dYk_ 

{2'kY 


' ( — 77 , 77 )^ 


f{k)b^^\k) 


d^k 


{2tt 


\d ■ 


(3.39) 


The functions Gz and D have these symmetries, so that we can replace e'^'* in (13.141) by 
b^^Yk). In Section 15.11 we show how to compute In,i{x), and in Section [5l2l we bound the 
other integrals in terms of In,i{x). 


Decomposition of the two-point function. We decompose Gz{k) and AGz{k) into sev¬ 
eral pieces, which we then bound in the next section. We start with some preparations for 
this decomposition. For the SRW-contributions, we define 

C*(k) = -.- - -.-. (3.40) 

l-F,(0)+a^,,[l-D(fc)] 
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As Fz{0) < 1 and ap^z > P p^ know that C*{k) < -^^—C{k) < (3 ^^C{k). Farther, we 
conclnde from 


C*{k) = 


1 


1 


“F,: 


1 ~ -^^( 0 ) + OiF,z 1 ~ ^ 3777 ^ 7 ^“- 1 “ 

1 


«F,; 


1 — Fz{^) + Of 


P'z ( 0 )+aF ,2 ^— P ' z { 0 )-\- o , 

Cx(^z) (^)? 


[1 - D{k)] 


with 


X{z) = 


1 


a 


F,Z 


2dl-Fz{0) + ap,z' 

and the monotonicity of Cx{x) in the parameter of the generating function A that 


C*ix) < 


1 — Fz{0) + Of,; 


Thus, C* can be bounded by ap^zC in x-space as well as in A:-space. 
We abbreviate 


-C{x) < - C{x) < — — C{x). 


a 


F,Z 




(3.41) 


(3.42) 


(3.43) 


Ff,z{ 0] k) 

m 


RF,ziO) ~ RF,zik), 

Rp,z{0-,k)C*{k) 
l-Fz(k) ’ 


R ‘ S >, z ( fi ] k) — i?#,2;(0) — R ^^ z { k ), 

M*{k) = b{k) - 2b^^^{k)C*{k). 


(3.44) 

(3.45) 


Bounds on key quantities. To bound the remainder terms, we define 

Kaf = z - \ -> so that --1— < K^pC{k). (3.46) 

!L,f-Par,f l-Fzik) 

Further, we assume that f 2 {z) < F 2 , i.e., that 

\a.{k)\ < ^^r2C(t). (3.47) 

and from now on abbreviate F 2 = We use the bound (|2.12ll in Assumption 12.71 to 

obtain 


|i?F,.(0; A:)| < [1 - b{k)] Y, \\x\\l\Rp,z{x)\ < [1 - ^(A:)]/3a«,f 

X 

Arguing as in (I3.12I1 - (I3.13I) we can show that 


4k Rp^ziS^] k)\ 


^ ||x|||i?i., 2 (x)e'^'^| < ( 5 ar,p 


We conclude the same bounds for A:), where /3 a«,f is replaced by /3ar,<s>- 

Combining these bounds with C{k) = 1/[1 — b{k)] we obtain 


\m\ < 


Rp,z{0-,k)C*{k) 


l-Fzik) 


^ /^Afl.Fi^AF^ (^) — 


Par,f^af 


Cik), 


and 


1 


\[R^Xk) - RFA0-,k)Gz{k)]C*{k)\ < - i/3p,^ + /3AR.Fb2)C{k), 

^ F.Z 


(3.48) 


(3.49) 


(3.50) 


(3.51) 
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Decomposition of /\Gz{k). We decompose AGz{k) into five contributions Hi{k). The 
dominant contribution is Hi{k), which is defined to be 

Hiik) = + a^^zb{k))C*{k) + a*,,) C*{k)M*{k). (3.52) 


The remainder terms H 2 {k)-,H^{k)^Hi{k) and H^{k) are defined as 


Mk) 

Hsik) 

Hi{k) 

H^{k) 


— — ( 0'F^ziC'i,z + Oiif,^zb{k)) (C*{k) + 


l-Fz{k) 


+ 


^z^E{k)M*{k) 




=2 


{l-Fz{k)f 

b^^^{k) 


{aF,zGzik) + iE{k) - 


1 - Ezik) 

AR ^^ z { k ) ARp ^ z { k ) 


OtF,Z ~ 1 

l-Ez{k)J ’ 


Gz{k), 


= - 2 


l-Ez{k) l-Fz{k) 

Tfs=iidsRFAk)? + 2ap,zdsb{k)dsRFAk) 


( 1 -F,(fc ))2 


Gzik) 


{l-FzikWfA 


^ (dsRg>,z{k)aF,zdsb{k) +ds^z{k)dsRf 


In Appendix [Cl we explicitly show that 


(3.53) 

(3.54) 

(3.55) 

(3.56) 


5 

-AGzik) = (3.57) 

i=l 

This computation is quite long and tedious. However, since it is crucial to our analysis, we 
give the derivation in detail in Appendix O Let us now give some insight into the origin of the 
different contributions. The first term Hi{k) is a SRW-like contribution that can be bounded 
similarly as in Section [3.3.31 The second term 7^2(^) corresponds to everything that has the 
factor M*{k) and a remainder term. In Ef^A), we collect the remaining Z)®™(A:) contributions. 
In Hi{k), we put the contributions of AR^Ak) and ARpAk), and in H^{k), we collect all 
products of single derivatives. 


3.3.5 Improvement of bounds for /a 

In this section we bound rA{x) by deriving bounds on 


R-iix) = 


(—7r,7r)'^ 


Hi{k)&{k)GAk)b^'^\k) 


Ak 


id’ 


(3.58) 


for z = 1,..., 5. We do this for z = 1, 2,... , 5 one by one, starting with 77”’^(x). This is 
the most technical part of the analysis. We bound each term of (x) using the bounds of 
Assumption 12.71 and the SRW-integrals (|3.35p - ()3.38p . 
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Step 1: Bound on We first recall the rearrangement of (I3.28p - (I3.30I) to see that, 

for m > 0, 


r (jdu 

/ D{k)^C*{kr+^M*{ky^-^-^^ = Y, \\y\\lC\y){D*^ * {cy^){x - y) 

J ( — 7 r, 7 r)^ y 

< \\y\\lC{y){D^^ * - y) 

y 

=(«.,.(3.59) 
We perform the bounds for n = 0, n = 1 and n = 2 separately. 


Bound on a weighted line (n = 0). Substituting the definition of we see that 

(|3.52l) leads to three terms. We bound the first and second term using (I3.59p . For the third 
term we can not use an x-space representation like in (I3.59p . To bound this term we repeat 
(|3.27p - (|3.3Up to see that 


D{kyC*{k)M*{k)e 


ik-x 


d^k 


^(—7r,7r)^ 

^ _ hj^*i (J*)ix) + ^ 


d' 

ms 1 11 .^ 

< -/i,i+i(x) + ^^^/ 2 ,/(x + 2 e,). 


In this way, we obtain 

, 0 ,/ 


(3.60) 


<l3,^^Jo,iix) + fd\^,^\Jo,i+iix) + ^^/i,i+i(x) + X]^ 2 ,«(a: + 2e,). (3.61) 

— n / F ^ TP I 


Bound on a weighted bubble (n = 1). To bound F}’^(x), we expand Gz{k) as follows: 


Gz{k) ={c^,z + a^,Mk))C*{k) + {R^,zik) - F^,,(0; k)Gzik) C*{k), 


so that 

U 


l;i(x) = / (c*,. + a^,zD{k))Hi{z)C*{k)D^-\k)&{k) 

J (—7r,7r)^ 


d^k 

(27r)‘^ 


+ 


' {—7T,7r)^ 


R^Ak) - RF,Mk)Gzik)) C*{k)Hi{z)b^^\k)&{k) 


d'^k 

{2ttY 


+ 


/?' 


■a,F 


with Tn^i as defined in (13.3711 . 


(3.62) 


(3.63) 


We bound the first line using (|3.59p and the second line using (|3.5ip to obtain 

|f};^(x)| <^“^(^,,*)^Ji,z(x) Ji,z+i(x) 

+ (/3|q.<i>i)^^^ ^77o,z+i(a:) + ^^^{f3^a,A‘^Ji,i+2{x) 

I3r,<s> + /^aa.fF^ 


-AcfTsAA + /3|^.4.|F3,z+i(x) + /3|c,4.|T2,z(x)), (3.64) 
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Bound on a weighted triangle (n = 2). We decompose G‘l{k) into two terms as 


GM" =C*{kf{c^,, + a^,,D{k)y 


(3.65) 


+ 


R^,z{k)-RF,z{^-^k)Gz{k) C*{k) {c^,, + a^^Mk))C*{k) + G,{k) . (3.66) 


We compute the contribution (|3.65p to be 

Hi{k)C*{kf{c^,, + a^,,D{k)? 

= + a^^zD{k))‘^ [c(f,z{c,s,,z + a-s>,zD{k))C*{k) + a^^z\M*{k)C*{k^. (3.67) 


We expand the brackets and then use (j3.59p to bound this contribution by 
^ Hi{k)C*[kf{c^^z + a^,zD{k)fD^^\k)&{k)- 


(3.68) 


^(—7r,7r)^ 

— \Pc,'i<^2,i{x) + /3\a,f\Ji,i{x) + 3/3|„,<i,|772,i+i(x)] 

+ + 'iJ2,l+2{x)\ + ^(/^C,*)^ [J2,l+3{x) + Jl^l+2{x)] , 

where we use that Jn,i{x) > 0 for every x € n, Z > 0 by p.28P " ()3.30l) . 

Using P3.51I) we bound the absolute value of the minor contributions given in (I3.66p by 


~k / (/^e,$ + P\a,‘S>\\Dik)\) p/ \ C [kS^ 

\ ^ F,Z 


a 


(3.69) 


F^Z 


Thus, we bound the contributions due to (I3.66P by 


< 


/ I77i 

J ( —TTjTt)^ 

( Pc,i 


F,Z 


OiF,Z 


+ l^\a, 


2 1/3 

\—a,F 

Pr,<s> T Par ,^^2 


P 

—a,F 


+ r'o 


Pc,<s>T4,,i{x) + /3|,,<s|r4,/+i(x) + /3 |„,s|T3,z(x) 


(27r)'^ 

(3.70) 




—a,F 


P^,f,T4^i+i{x) + /3|„,s|T4,z+2(a;) + /3|c.4.|T3,«+i(a:) 


Conclusion of Step 1. We have bounded the contribution due to Hi{k) and have obtained 
that 




'(IMO 
< (TOIl 
(lOTll + (13^ 


for n = 0, 
for n = 1, 
for n = 2. 


(3.71) 


By the sum of two equation numbers we mean the sum of the terms given in the right- 
hand sides of the corresponding equations. As for z = zj, this bound uses 1^+2,px) and can 
therefore not be used in d = dc + 1, dc + 2. We have chosen to use these bounds, even if other 
bounds would be available, as they give numerically better bounds. 
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Step 2: Bound on For this bound we use 


Kii^) = 


' ( —7r,7r)^ 


\&{k)\C{k)^\D^^\k)\\M*{k)\ 


df^k 

W) 


d ’ 


(3.72) 


which is an adaptation of defined in (|3.37l) . In Section [Ql we will bound T*^ in the same 
way as T^^i- We bound the absolute value of H 2 {k), defined in (I3.53j) . by 


+ /3|c.,*i )Par,fILaf ^|iV7*(A;)| 


\H 2 {k)\ <[a^A /3,,^ + /3K*||^(fe)|)(a-i +K, 

+ ^^,FPAn,^KlFC{kf\M*{k)\. 

We use that \Gz{k)\ < T 2 C{k) and the integrals T*i to bound by 


Of 


(3.73) 


|772:.(^)I < /3A«..KA,(r'2)"[(/3„,T*+2/x) + /3|.,,|r„V2,z+i(x)) +^Ar-j (3.74) 

+ /3|„,,l^;i^r*+i,,(x)l +P^,,/3An,.Kl,T*^2Ax). 


Step 3: Bound on 7^3 2 (x). In ()3.54p . we have defined H^^k) to be 

Hs{k) =2-^^^ (a^,.Gz{k) + a*,,) (E{k) - ^ 


l-Fz(k) 


l-Fzik)J • 


We bound |i73(A:)| as 


i73(fc)| < 2D^'^{k)K^^C{k) (a^,zT'2C{k) + a^A K^FC{k) 


f^AR,I 

OiF,Z 


and use this bound and the integral Un,i, defined in (I3.38p . to bound |773 2 (x)| as follows: 


(3.75) 


+ |a^,,-I| . (3.76) 


|77"j(x)| < 2(r')"+^^L {Par,f + max{|/3„,, _ i|, ^ _ i|}) ( 3 . 77 ) 


\n+l 1^2 


•a,F 


+ 2 (r'+ niax{|/3„,^ - 1 |,\8 - 1 |} Un+ 2 A^) 


Step 4: Bound on TCI’1{x). We first bound HAk) in Fourier space 


as 


\HAk)\<KAF (/3 ah,<. +/ 3 ah,^F '2 


Then, we use the definition of in (|3.36p to bound 

\KiA)\ <Kaf {f 3 AR,^KnAA + f 3 AR,Fr' 2 Kn+iAx)) 


(3.78) 
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Step 5: Bound on T-Q'l{x). We recall that 


Mk) = - 2 


Eti dsRpAk){‘^ap,zdsD{k) + dsR^Ak)) 
{l-F,{k)f 


{l-F,{k)Yjz 


G,{k) 


^ [^dsR<t,z{k)aF,zdsD{k) +ds^z{k)dsRf 


(3.79) 


S=1 


To bound the single derivatives we note that for a totally rotationally symmetric function /, 
see Definition 12.51 the following holds: 


dsf{k) = \ '^Xsf{x)F^'^ 
X 

for s G {1,..., d}, so that 


f{x)xs sm{ksXs) cos(/cp 

X l/yts 


(3.80) 


\dsf{k)\ < ^ |/(x)||x^sin(A;5Xs)|. 
X 

Since | sin(nt)| < n\ sin(t)| for integer n, we obtain that 

\dsf{k)\ < |sin(A:^)i^|/(x)|x^ 

X 

The total rotational symmetry of / also implies that 


(3.81) 


(3.82) 


x)|||x||2 


(3.83) 


for all s, t G {1,... , d}. From this we conclude for two totally rotationally symmetric functions 
/, g that 


\8j{k)d.g{k)\ < Eti Er Ey |9(!/)|!/," 


S=1 


= ^"“WEJI^IIil/(^)IEJ|y|lil5(y)l, 


(3.84) 


where we recall (j3.23p . Using this relation we can bound H^{k) by 


\Mk)G^z{k)\ < 2Kl^T'^+^C{kr+^b^'^{k){2aF,zP^u,F + Pln.^ 

+ ‘^I^^F^2^G{kb~^‘^b^^^{k){aF,z/3AR,i' + l3\a,^\l3AR,F T /^Afl.F/^Afl,*)- (3.85) 

and obtain the following bound on 1-Q'\{x)-. 


|77^;'(x)| < 2KLr'’^+i(2/3„,^/3^H,^ + /3E,^)C/„+3,z(x) 

Afl,« + ld\a,ii\(d AR,F + /^Afl.F/^Afl,* )Un+ 2 ,l{x)- 


(3.86) 
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Final bound on f^. In this section, we have bounded /s by 


fsiz) < max 
{n,l,S}€S 


sup^gs I (l3TTi) + (ICT) + (I37m + (137781) + (13:861) | 


Cn,l,S 


(3.87) 


We recall that by the sum of several equation numbers we mean the sum of the terms given 
in the right-hand sides of the corresponding equations. 

In summary and recalling Definition 12.91 when P( 7 , r,^) holds, this bound on is 

smaller than 73 . Thus, the improvement of all bounds is successful and we have thus success¬ 
fully performed the bootstrap. The computation of a numerical value for the bound in (I3.87p 
requires the computation of the SRW-integrals In,i, Kn^i,Tn^i,Un,i- In Section IT^ we show 
how to bound SRW-integrals and explain for which x the supremum over S is obtained. 

The bootstrap function provides various bounds on weighed diagrams. The real size 
of these diagram depends heavily on the values of n, I and the set S involved. For example, 
we can expect that 'H^z^{x) is of order 0(1) while 'H^z^{x) is of order 0{d~‘^). Since the form 
of the bounds on 'H^'\x) is the same for all n,l, we have introduced the constants Cnj^s to 
merge them into one bootstrap function. Alternatively, we could consider to consist of 
multiple bootstrap functions that are individually bounded by within the bootstrap 

argument. 


4 Rewrite of the NoBLE equation 

In the preceding part of this paper, we have performed the analysis using the form (|1.37p for 
the two-point function. This form is related to the classical lace expansion. We have decided 
to use this as it considerably simplifies the presentation of the analysis in the preceding 
section. 

In this section, we first derive this characterization from the NoBLE equation, meaning 
that we identify oi^^z-,oif,z-,Rf,z and R^^z- Then, we translate the assumptions made on the 
rewrite (I1.37P into assumptions on the NoBLE-coefficients Hi,, 

The aim of the rewrite is to extract the dominant SRW-like contributions from and 
Fz, see ()1.33p - (|1.34l) . These SRW-like contributions will give rise to 0 :^, 2 , aF, 2 , £*, 2 , The 
remainder is put into i?F ,2 and R^^z- 

Here we show how we extract SRW contributions from ^z and Fz and use them in our 
analysis. More terms could be extracted from and Fz, thereby reducing the value of i?F ,2 
and R^^z and thus increasing the performance of the perturbative technique. This might allow 
to prove the infrared bound in even smaller dimensions above the upper critical dimension. 
We however found the possible gain not in relation with the necessary efforts. 


4.1 Derivation of the rewrite 

In this section, we rewrite the functions ^z{k) and Fz{k), as defined in (ll.33p - (ll.34l) . and 
identify a.j,^z,ctF,z, Rf,z and Rs, 2 - The NoBLE-coefficients are defined as alternating series of 
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(4.1) 


non-negative real-valued functions 

oo oo 

2.W = j;(-l)'"2f>{x). S,{x) = j;(-l)^2f>-‘(,T). 

7V=0 iV=0 

oo oo 

nix) = nr(x) = j;(-i)^nf>"(,T). 


(4,2) 


N=0 


N=0 


4.1.1 The model-dependent split of the coefficients 

When rewriting the two-point function, we extract a major SRW-like contribution. We are 
guided by the intuition that coefficients are of order 0{{2d)~^) and that the main contributions 
to the NoBLE coefficients are 


=?>(ei) =» =«(ei) (4.3) 

=<">-‘(e.) f..nm'‘>'‘(e.). (4.4) 


Due to the limitation of our bounds it is not beneficial to extract all these contributions. 
Thus, we create a model-dependent split of the coefficients to improve the performance of the 
technique. We define non-negative functions 


'='(0) 

'='(1) 

^ aA,z'> 

^ Q.,II,Z 5 

^a,7,2? 

^ (y,ii,z ? 



tt(0),z,,«: 

^^Oi,Z 




^OlAI^ZI 

^(0) 

^(1) 

^ R,I^Z^ 

iTrCO),/- 

^ R,II,Z1 

^ RJ^Zi 

^ RJI^ZI 



^^R,Z 


^R,Z1 

^R,I,Z1 



Here these functions satisfy that, for = 0,1 and all x € 






(a;) = + ^AiiAA, 


X = 


■^°^’'z(x) -h eAAA = 


-‘aA,z \ 


-‘a.ii.z 


A) + 


and, for x G Z'^ with llxllo > 1, 


■'^^(x) = 0, ^aAzA + d) = 0, ^aAAx) = 0, 


AAzA + d) = 0 , 


,(N),t 


(N),t 


Hx) = 0 , 


and 


uA/Ax) = 0 , 


for X ^ {ct, Ct -)- Ck}. Further, these functions have the same symmetries as the original coef¬ 
ficients. The idea behind these two different splits (giving rise to the terms with subscripts 
I and II, respectively) is that we split off specific contributions that can be explicitly incor¬ 
porated in the constant and D{k) terms in our expansion. Contributions with subscript I 
correspond to x for which ||x —ej| < 1, while contributions with subscript II correspond to x 
for which ||x|| < 1. In Fourier space, this corresponds to contributions with a factor 
and e^^^, respectively. See ()4.14p below for how such contributions will arise. 
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4.1.2 The Fourier inverse of F and 4> 


Throughout this section, we omit z from our notation and write, e.g., = /x and Fz{k) = 

F{k). As a first step, we use the Neumann-series to rewrite F and $ into a form without 
matrices. We use that (D(A:) + = (D(—A;) — /xJ)/(l — to rearrange F{k) as 


F{k) = 


/X + 

1 — /x^ 


1 + 
1 + 


B{k) + fi3 + n{k) 

1 


1 -1 


1 + 


1 — /x^ 


(D(-A;) -;uJ)n(/s) 


(D(-A:) -/xJ)l 




n =0 


1 


E E 1 + 


I — fF 

(-ly 


(D(-A;) - fiJ)n{k) {t){-k) - /xJ)l 


1 — ^ ^ V V (1 — FF 

^ n=Oio,...,in ' 


(4.5) 




We define F^ as the nth contribution in the sum in (14.51) and analyze these terms separately. 
The Fourier inverse of Fn = Fn^z is given by 


Fo{x) = ^ ^ ^ (4,-e, - ^JFo,x + + ej - /x4''(x)) 


Fn{x) =^^ Y 

tOi-.-i'-n Xi'FZi 


(4.6) 


(-l)"(4o,o + ^^°(xo)) 
(1 - /x 2 )«+l 


^n—1 


- /xn-‘-i’^«(x,)) 


^s=l 


X + e,J - /xn‘"-i’^-(x„ + 

- /xn-''"-i’^"(x„ + e,J + ^2n-''"-i’^"(xn)). 


(4.7) 


In a similar way, we define such that 


^{k) = Y^ ^nik), so that also <I)(x) = Y^ ‘l’n(a^)- 


n =0 


n =0 


(4.8) 


These function are given by 


ioW =1 + m - E(‘ + 'i'‘(»:))(^We-''=- - f=-‘(t)), 

i — l 

E (i + 'i'nt)) .), Wi 

L 0 ,...,Ln ^ ^ ^ 

n 

X (^n^^-i’‘-(A;)e“'*^‘-i -/xn-^“-i'^“(x,)) (H‘"(A;)e-'*^‘- - fiE-^^{k)). 

S = 1 ^ 


(4.9) 

(4.10) 
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The Fourier inverses of these functions are 


^>o(x) =(5o,x + H(x) - ^ ^ y'(^o,y + -y + e,)- nE ‘'{x - y)), (4.11) 

i-,y 

<^>nix)=y Y (<^0,xo + ^'°(3:o)) _ ^2\n+l 

'-Ovi^n Xi:J2i Xi=x 


S =1 


4.1.3 Definition of the rewrite 

In the rewrite, we extract explicit terms that are independent of k and terms that involve 
D{k) for F and $. Everything else is put into the remainder terms i?F,z and R<s,^z- The major 
contributions that we can extract are part of Fq and $o- Also Fi gives some contributions. 
We begin with Fq and rewrite it as 

Mk) + *'('=)) («■“■ - ") 

^ i 

OO 

=-^(2dD{k)-2dy+Y -^)). (4.13) 

^ Af=0 t 

Recall that the lace-expansion coefficients are defined via an alternating series of non-negative 
functions, see dEI, (gJl). For N = 0,1, we split the sum into 

^ T(^)’‘(A:)(e-‘^^ - y) = Y + ej - /iT(^)’^(x)) 

i L^X 

+ 2db[k) Y {^ax\ei + e.) - 

K 

+ YY e'"'" + ej - h'I^S’/(a^)) , (4.14) 

k xGZ'^ 

where we see how the splits involving the subscripts I and II are used to extract random-walk 
contributions. From 

= ( 1 ^ 2)2 E (l + ^'°(fc)) (4.15) 

we extract the contribution of 1 x x e“‘^‘i and split it as 

^e-i(fc.o+fcu)n(“)’‘0’^i(A:) = 2dY (n((’^’^’''(ei + e«) + i)(fe)n((’^4A(g^)^ 

LQ^il K 

+ YY + e,J. (4.16) 

'■O,'.! a;gZ‘* 

We have now collected all the terms of the split in the lines (I4.13I) - (I4.16I) . The constant terms 
contribute to Cp^z- Terms involving D[k) give rise to ap^z- All other terms contribute to 
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Rp^zik). Thus, we conclude that 


Cp.z — 


2(1^1 2dfi 


+ 


1- ni 1 - Ai 




AfG{0,l} 


2d^z 


Q^F,; 


2dfj, 




(4.17) 


i-m^l 


Afe{o,i} 


2dH:z 




^nsw.’‘(ei), 


and 


(4.18) 

RF,z{k) = F^{k) - Cp^z - aF^zi){k), (4.19) 

which is the sum of the final contributions to the right hand side of (I4.13p . (|4.14p . p4.16p and 
the remainder of (I4.15p . We rewrite 4>(fe) in the same way. We begin by noting that 


Mk) = l+ E (-1)'^HW(A:) +J](-l)^H(-)(fc) 

AfG{0,l} 


(4.20) 


N=2 






1 — /i^ 

For = 0,1, we split E^^\k) as 

E^^\k) = HW(0) + 2dDik)E<~^\ei) + E^^\k). (4.21) 

Further, we extract the contribution of the factor 1 and in the second line of (|4.2UI) as 


^ ei^-(H™’‘(x + ej - ArH(°)’-X0)) =2d{Eff{e,) - mH™;)(0)) 


(4.22) 


+ 2db{k) ^ {E^°]t{en + e,) - /uH™i)(eK)) 

K 

+ E E + e,) - ^^Ef:-\x)). 


We define 


iVG{0,l} 






( 0 )), 




2d Y. (-i)''2S(ei)--i%E{=«(*»+'■) 

Afe{o,i} K 


/x^“™;)(eK)) , 


Ri>,z{k) =^z{k) - c^,z - a^,zD{k). 


(4.23) 

(4.24) 

(4.25) 


This completes the derivation of the rewrite (11.351) and pi.361) and identifies aF,z,o:.i>^z,Cp^z:, 
c.i>^ 2 :? Rp^z nnd Rqt^z' 

At this point, it is worth mentioning that this is not the only possible split. Indeed, we 
could try to put more terms into aF, 2 ; 0 «, 2 , thus reducing Rf^z,R'S>,z, and thereby improving 
the efficiency of the analysis. However, numerically we found that the possible gain would 
not be in relation to the necessary efforts, so we refrain from this. 
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4.2 Assumption on the NoBLE coefficients 

In this section, we reformulate Assumption 12.6112^81 on 0^,21 ct «,25 Rf,z and in terms of the 
NoBLE coefficients. We assume that the NoBLE coefficients have the following properties: 

Assumption 4.1 (Symmetry of the models). Let t, k € {±1, ±2,... , ±4}. The following 
symmetries hold for all x € z < Zc, N ^ N and i, k: 

For all A G N, the coefficients 


E=?’'‘w Enr^w, (4.26) 



L 

L 

L,K 


as well as 

the remainder terms of the split 





(x), Y 4'k"!n,2 {x), 

L L 

Y 

L 

L 

Y^'A^’Ax), 

L^K, 

(4.27) 


are totally rotationally symmetric functions of x ^ Z*^. Further, the dimensions are exchange¬ 
able, i.e., for all l, k, 

i5f)-^(0)= = (4.28) 

k' l' 

The next assumption states a bound on and in terms of 

Assumption 4.2 (Relation between coefficients). For all x G Z'^, p < pc, A G N and 
t, K G {±1, ±2,..., ±d}, the following bounds hold: 

XI,(N)A(^) <^Hf>(x), nf>-^’'^(x) < (4.29) 

As explained in Section IXTl to successfully apply the bootstrap argument, we assume that 
the coefficients obey certain bounds when the bootstrap assumption fi{z) < Ej holds for all i G 
{1,2, 3} for a given z G [0, Zc)- These bounds do not depend on the value of z. However, their 
form is delicate and depends sensitively on the precise model under consideration. We assume 
that the same bounds hold for zj regardless of the values fi{zi), f 2 {zi), f 3 {zj). Assumption 
14.31 is the most technical assumption of this paper, and is phrased so as to allow maximal 
flexibility in the application of the NoBLE: 

Assumption 4.3 (Diagrammatic bounds). Let ri,r 2 ,r 3 > 0. Assume that z G {zi,zf) is 
such that fi{z) < R fori G {1,2,3} holds. Then Gz{k) > 0 for all k G (— 7 r, 7 r)'^. There exists 
— 1)/^ >0 such that 

— /i 2 > /3 . (4.30) 
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Further, there exist o>/^a2‘ — 0; such that 


e(N)(o) < ;3 (n)^ 
||x||i=<">(i) < A'2, 

X 

Y1 11^ “ 


y: iixiiiHr’^(x) < 


for all N > 0 and k G (—tt, vr)'^. Moreover, we assume that Yl'N=of^»^^ ^ °° 
{H, AS, {AH'', 0}, {AH'', i}} and that 


{2d - l)fl, 

1 hz 


Es 

Ar=0 


(N) 


< 1 . 


Further, there exist such that 

’ £ 1 ® ’ —En 

^(o),.(o) > /5(o)^ 


^n(i),.A(o) > . 


Additionally, there exist ^ with 

<S™(0)-HW(0)</3^:;>^, 

and jB^S rn^ 0 suc/i that 

K 

K. 

Also, there exist ^^ehc’ 

~'^E®L^ — ^ “ '^amziei + Ck)) < 

K 


^Eo. 


(0) 

E Hq 


For N = 0,1, there exist P^Xi, /^aXr,ii - 


?(N) 


(N) 


3(N) 


?(N) 


E^Ew <«. 

X 

Ell* 

X 

WlE^lix) < X 

Y'^wXzix) </3i%, 

X 

Ell 

X 



E 

lklli^kE,^(2;) 


:,7?’ 


, 7 ?, 7 ’ 


(4.31) 

(4.32) 

(4.33) 
for • G 

(4.34) 

(4.35) 

(4.36) 

(4.37) 

(4.38) 

(4.39) 

(4.40) 

(4.41) 

(4.42) 

(4.43) 

(4.44) 

(4.45) 
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Further, there exist 


I^As\r,ii’ /3n%; /^An.fl - 


X 


Y 11 ^ “ -h ej < /JaLgh./j 

X 

(4.46) 

X 

</o(0) 

/ ^ — PAs^rjV 

X 

(4.47) 

X.L 

<q(0) 

Y Il^ll 2 n™/’''(a; + e, + e«,) < (IYr- 

X.L.K 

(4.48) 


For all • G {H, S'", AH, {AH'", 0}, {AH'", t}} and A G N, depends only on Fi, r 2 , Ta, d and 
on the model. If Assumption \2.S\ holds, then the bounds stated above also holds for z = zj 
with the constants /3, only depending on the dimension d and the model. 


Only the bounds (|4.30p - (|4.33p are essential to perform the analysis for the NoBLE. The 
bounds stated in (j4.36l) - (l4.48p are used to obtain good bounds on c,s,^z,0(.j>^z,CiF,z, Rf,z and 
that allow us to increase the performance of the analysis and to show mean-field result 
in lower dimensions than otherwise possible. 

We denote by and Z?®™" the sum over all (resp. odd/even) N of (li^\ i.e., 


CXD 

/3r^ 


oo 

^odd ^ ^ ^(2N+1) 

N=0 


^even 


CXD 


E -Si”’. 

N=0 


(4.49) 


for • G {H, AH,{AH^0},{AH^i}}. By (gT]), the values and (-,5°?'^) 

are explicit upper and lower bounds on H^(0) and E(,(0), respectively. By Assumption 14.21 
they also imply bounds on T^(0) and n^'^(O). 

We next discuss the left-continuity of the coefficients at z = Zc- 

Assumption 4.4 (Growth at the critical point). The functions z i—>■ 'Ez{k),z i—>■ H‘(/c),z i—>■ 
^f{k),z Ilz'^{k) are continuous for z G {0,Zc). Further, let ri,r 2 ,r 3 > 0 6e such that 
fiiz) < Fj and that Assumption \4.d\ holds. Then, the functions stated above are left-continuous 
in Zc with a finite limit z ^ Zc for all x . Further, for technical reasons, we assume that 

Zc < 1 / 2 . 


In the remainder of this section we show that Assumptions 14.1114.41 imply Assumptions 
12.6112.81 as formulated in the following proposition: 

Proposition 4.5 (Translation of the assumptions). The assumptions stated in Section\2f^ 
are implied by the assumptions stated in Section \^T^ More precisely, 

(i) Assumption \4.1\ implies Assvmvtion \2. (A 

(ii) Assumptions \4.1\ - \4-^ imply Assumption \2.7[ 

(in) Assumptions 14-21 - 14.41 imply Assumvtion \2.8l 

The proof of parts (i),(iii) are relatively straightforward and are performed in Section [4.31 
Part (ii) is proven using a tedious, but also straightforward, application of the bounds stated 
in Assumption 14.31 We add the details of this in Appendix [Dj 
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4.3 Properties of the rewrite 

In this section, we prove Proposition I4.5l ii and (hi). 

Proof of Proposition \4-5^ i)- In Assumption I4.1( we assume that the dimensions are inter¬ 
changeable (recall (|4.28p l. so that (j2.7p clearly holds. Further, we also assume in Assumption 
o that the two-point function x i—is totally rotationally symmetric. 

To see that Rp^z and are totally rotationally symmetric, we note that convolutions 
maintain symmetry and that the NoBLE-coeflicients, in the way they are combined in the 
definition of and R^^z, are thus also totally rotationally symmetric. This completes the 
proof of Proposition I4.5l il. □ 


Proof of Proposition \4-5^ iii)- We prove the statement in three steps: 

(a) We prove that [D(A:) -|- is well defined for all k and z < Zc, 

(b) We conclude that z ^z(k), z Fz{k) are continuous at z and well defined at Zc, which 
implies that the bounds stated in Assumption 12.71 also hold for z = Zc] 

(c) We show that Gz{k) can be continuously extended to z = Zc for 7 ^ 0. 

This implies the desired statement that z Gz{k) is left-continuous at z = Zc for any A: 7 ^ 0 
and that the bounds stated in Assumption 12.71 also hold for z = Zc- 
(a) We begin by showing that 


[D(A:)+;U,J] ^tlz{k) 


sup max 

lklloo=i " 


([D(fc) + /izJ] ^tlz{k)v 


< 1 . 


(4.50) 


We start by noting that [D(fc) -|- /x^j] ^ ^ (D(—A:) — ftzJ), so that 


([i)ik) + f,zJ] "Uzik)v) (4.51) 


Thus, for V with ||u||oo = 1, 


[-b{k) + fizJ] "flz{k)v 


^ 1 + Mz II I 

- rr?"”' 






< 


1 - Aiz 


N^K,,X 

~ 1 , 


(4.52) 


1 -1 


which proves (14.5011 . From (14.501) . it follows that the matrix I -|- D(A:) -|- /x^J n 2 (A:) is 
invertible. Then, we use standard linear algebra to compute 


I + [D(A:) +/XzJ] ^flzik) ‘ [D(A:) + A^zJ] ' [D(A:) + ^ + n,(A:)] (4.53) 


-1 ^ , 


-1 


I + [D(A;) +/XzJ] ^nz(A:) ‘ I + [D(/j) +/x,j] ^II,(A:) 


-1, 


= 1 , 


which implies that the matrix D(A;) -|- AizJ + nz(^) is invertible. 

(b) By Assumption l4.4l we know that the NoBLE coefficients are continuous in z. This also 
implies that D(A:) -|- a^z J + nz(^) is continuous in z and, as it is well defined by (a), its inverse 
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is also continuous. Further, we note that the bounds on the coefficients /?, are independent 
of the value of 2 ; € {zj,Zc) and the coefficients are left-continuous in 2; = Zc- Reviewing the 
definition of ^z{k) and Fz{k) in ()l.,'13p - ()l.,'14j) . we conclude that these functions are continuous 
in 2; G [zi, Zc) and left-continuous at 2; = Zc- 

(c) The dominated convergence theorem implies that 


Gz^ix) = lim Gzix) = lim f 

z/'ZcJ(. 

Kik) 


^zik) 


Ak-x 


i 


dA 


I - Fz{k.) (2vr)'= 
d<^k 


ik-x 


{-Tr,wy I - Fz^{k) (2vr) 


d’ 


(4.54) 


where we use the left-continuity of 2 i-)- ^z{k) and 2 i-)- Fz{k) at 2 = 2c proved above, and 
we further note that /3 ^ p ^ ^ ^ ~ 1 “ Pr,<s> > 0 together with (I1.37P 

imply the infrared bound holds uniformly in [ 2 /, 2 c) for $ 2 (A;)/[I — Fz{k)], so that the integral 
in (I4.54P is well dehned. Thus, we still have the Fourier representation (11.321) . with the 
understanding that Gz^{k) for k ^ 0 is dehned by Gz^{k) = $ 2 ^(/c)/(l — Fz^{k)). Since 
^2^(0) = 1, this characterization can not be used for k = 0. This completes the proof of 
Proposition I4.5f iiil. □ 


5 Numerical bounds 

In this section we discuss the ideas underlying the numerical computation of our bounds on 
the NoBLE coefficients. These ideas are model independent, while the implementation itself 
is not. We hrst explain how we compute the numerical bounds on the SRW-integrals that 
we have used for the improvement of bounds in Section [3] and to obtain numerical bounds 
on the coefficients. Then, we explain how the bootstrap functions are used to bound simple 
diagrams. At the end of this section, we explain how we compute the as sums over 
which is not straightforward as the bounds on the NoBLE coefficients are stated in the form 
of matrix-products. 

5.1 Simple random walk integrals 

We bound the SRW-integrals In,i, Kn,i,Tn^i,Un,i dehned in (I3.35p - (I3.38I) . We hrst compute 
In,mA) then show that the other integrals can be bounded in terms of it. We compute 
using 


In,mix) — I-ii^m—lix) In—l,m—lix), (^-f) 

which is obtained by writing Dik) = 1 — [1 — Dik)] in (I3.29p . Using (15.11) . the problem of 
computing In,m for general n, m G N simplihes to the computation of In,o and /o,m for all 
n, m G N. 

5.1.1 Computation of the Green’s function 

We compute In,o in the same way as Kara and Slade in m Appendix B], as we explain now. 
Let bin, s) be the modihed Bessel function of the hrst kind and F(t, d, n) the modihed Bessel 
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function, i.e., 


OO (-»; I 

^ / c\2k-\-n 

i.(n,s) = 


fc =0 


k\T{n + A: + 1)! 


F{t,d,n) = e (5-2) 


see e.g. [19l (8.401) and (8.406)] or [TJ Section 9.6]. Using 

1 1 


[l-D{k)]'^ (n-l)!7o 


/*oo 

Jo 


(5.3) 


we compute 




1 ^ 

(^Tziyr 

M ~1 


(5.4) 


see [271 Appendix B]. Most mathematical software packages, such as Mathematica, Matlab, 
and R, come with a method to compute the modified Bessel Integral. We have used Mathe¬ 
matica which allows to control the precision of the computation. With the built-in function 
we compute in d > 15 and hflix) in d > 18 up to a precision of 10“^°. To be able to 

compute these basic SRW-integrals in lower dimensions, we implement the algorithm given 
in 1271 Appendix B], where also a rigorous bound on the error is proven. This algorithm is 
based on a Taylor approximation of the Bessel function. 


5.1.2 Computation of the random walk transition probability 

The computation of lo^mix) is a purely combinatorial problem as (2d)"*/o,m(a^) = Pm{x), 
where Pm{x) is the number of m-step SRWs with wq = 0,a;m = x. The value of Pn{x) can 
be obtained by simple combinatorial means. As an example, we explain the computation of 
P 6 ( 0 ). 

When cjo = wg = 0, the walk uses at most three different dimensions as it needs to undo 
all its steps. In the following, we distinguish between the number of dimensions used by the 
walker: 

O When the walk only uses one dimension, it steps three times to the positive direction 
(right) and three times to the negative direction (left). As any combination of left and 
right steps is allowed there are 6!/(3!3!) different possibilities for that. As there are d 
choices for the dimension used, there are d^ such walks. 

O When the walk uses two dimensions, it makes 4 steps in one dimension and 2 in the other. 
As any combination of moves is allowed, there are 6!/(2!2!l!l!) different possibilities for 
that. Further, there are d choices for the dimension in which to take 4 steps, and likewise 
d — 1 choices for the dimension where 2 steps are made. Thus, there are d(d — 1)^ 
SRW 6-step loops using steps in exactly two dimensions. 

O When the walk uses three dimensions, then there are 2 steps in each dimension. There 
are 6! different orders for these 6 steps (including the back and forth steps in each of 
the three dimensions). Further, we have to choose 3 out of the d dimensions (without 
repetition). This gives a factor 5 i_ 
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This means that 


Pg{ 0) = d 


3,3 


+ d{d - 1) 


6 

2,2,1,1 


+ 


d{d-l){d-2) 

3! 


6 

1 , 1 , 1 , 1 , 1,1 


where the multinomial coefficient is defined as 



ml 

ki\k2l... ,krl 


(5.5) 


(5.6) 


For our analysis, we use the values of Pn{x) for n G {0,..., 20} for about 24 different values 
of X. We have implemented a program for this, and the algorithm can be found in the 
accompanying Mathematica notebooks (see also Section El where the Mathematica notebooks 
are described in more detail). 


5.2 Bounds on related SRW-integrals 

In this section, we show how to bound the integrals defined in (j3.36l) - (l3.38p . This section is 
an adaption of m Appendix B.l] by Kara and Slade, who computed numerical bounds on 
these integrals to prove mean-field behaviour for nearest-neighbour SAW in d > 5. 


Bound in terms of Vn We first show how we bound the integrals defining Kn^i{x), 

Un,i{x) and Tn,i in terms of In,i, as well as the related integrals Ln{x) and Vn,i defined by 


Ln{x') — 


' (—7r,7r)^ 


C{kTlD^^\k) 


2 


and 


Vn,l = 


r D{ky[D^^^{k)Y d^k 


(5.7) 


(5.8) 


We use the Cauchy-Schwarz inequality to bound Kni{x) and Uni{x) defined in (13.361) and 
mB by 


<[In,2mLn{x)y^\ 


Un,l{x) < [Vn,2lLn{xp^. 


(5.9) 


To bound defined in (I3.37p . we use ()3.32p and |d)®'°(A:)| < 1/d, respectively, to compute 


\M{k)\ < \D{k)\ + 2\&^{k)\C{k) < \D{k)\ + min 1^, ^C(A:)| . (5.10) 

This leads to 

< [ \&{k)\C{kr\M{k)\\D(^\k)\^ 

< Kn,l+l{x) +min^Kn,l{x),‘^Kn+l,l{x)^ . (5.11) 

Next, we discuss improvements for the bounds on Kn,i{x) and Un,i- As \D^^^{k)\ < 1/d (recall 
()3.32p i. we know that 


Un,lix) <^Kn,l{x). 


(5.12) 
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(5.13) 


For X = 0 and even I, we use D''°'>{k) = 1, > 0 and (|3.25l) to compute 

Un,l{0) = — (/n,z(0) — In,z(2ei)) . 

For X = 0 we can use a better bound for Kn^i in the form 

K in')if/is even, 

\< 4,/-i(0)^/^/n,z+i(0)^/^ if / is odd. 

Moreover, we use a different bound for / = 0. We note that, for n > 1, 

1 _ 1 I D{k) ^ b{kf 

[i-b{k)]'^ [i-b{k)]^-^ [i-b{k)\^-^ 

which implies that Kn,i{^) < Kn-i,i{x) + iFn-i,z+i(a;) + Kn,i+ 2 ix), and thus 

Kn,o{x) <Kn-l,o{x) + [4-1,2(0)F„_i(x)]^/ 2 ^ (0)L„(x)]^ 2 , 


(5.14) 


(5.15) 


(5.16) 


Computation of Ln- By the definition of b^^\k) in (|3.34|1 . it is not difficult to see that 

{b^^\k)f (bk 1 


Ln (x) — 


[1 - b{k)]- {27T)d 2dd\ 


E E In,o{x+p{x]V, (5)). 


(5.17) 




The set Vd and the operator p(x; z^, 6) are defined in Definition 12.51 As we can compute 
In,o{x), we can also compute the sum in (j5.17l) directly. 

The value of In,o{x) only depends on the number of entries that have a given absolute 
value, so that we can reduce the domain over which we sum. We explain this in two examples: 


Example 1: Computation of L„(ei). As the first example, we show that 

Tn(ei) - d —(5.18) 

By symmetry. In,o(ei +p(ei; u, 6)) = /n,o(ei + 62 ) for all <5 G {—1, l}'^ and u ^Vd with ui / 1. 
This explains the third summand of ()5.18p . where we note that there are (d — l)!(d — 1) 
permutations u with vi ^ 1 . That leaves (d — 1)! permutations v with z^i = 1 . As all entries 
of p(ei; V, 6) except the first one are zero, the values ( 52 ,..., do not affect the summand. If 
(5i = 1, then ei +p(ei;z^,(5) = 2ei and if (5i = —1, then ei + p{ei]i',6) = 0. The two cases 
correspond to the first and second term in (|5.18p and complete the proof of (|5.18p . 

Example 2: Computation of T„(ei + 62 ). As the second example, we derive that 

^ , , (d-2)(d-3)^ , 

Tn(ei + 62) ——- 7 „^o(ei + 62 + 63 + 64) 

+ 2 d[d — 1) -/^ra,o(2ei + 62 + 63)) 

+ ^ (-^n,o(0) + 7n,o(2ei + 262 ) + 2/n^o(2ei)). (5.19) 
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There are 2{d — 2)\ permutations u with {^ 1 ,^ 2 } = {1)2}. Further, there are 2{d — 2)\{2d — 2) 
permutation u that map 1 to {3,..., d} and 2 to (1, 2}. That leaves 

d\ -2{d-2)\- 4(d - 2)!(2d -2) = {d- 2)\{d -2){d- 3) (5.20) 

permutations v that do not map 1 and 2 to the first coordinates, i.e., 12 for which { 1 ^ 1 , 2 ^ 2 } Fl 
{1,2} = 0. For these u, 

In,oiei + 62 + p(ei; i/, 5)) = /n,o(ei + 62 + 63 + 64 ), (5.21) 

which yields the hrst summand of (I5.19p . The second corresponds to the case that either ei 
or 62 is mapped to one of the first two coordinates. For example, let us assume = 1 and 
1^2 = 3, then 


61 + 62 + p( 6 i; V, J) € {2ei ± 62 ± 63 , ±62 ± 63 }, (5.22) 

depending on the sign of (5i. This gives the second summand of (I5.19p . If we map 61 + 62 to 
both of the hrst two coordinates, then 

61 + 62 +p{ei]iy,6) € {0,261,262,261 + 262 }, (5.23) 

which only depends on Ji and < 52 . This gives the third summand of ()5.19p and completes the 
proof of (|5.19p . 


Computation of Vn,i- The eqnality ()3.25p implies that 


D«“(A;)2 = ^^[1 - D{2k)f = ^ 


{2dy 


{2dY 

From this, we conclude that 

b{ky[b^'^{k)f d<^k 




,j2i(A:i+A:„) 


{2df 


(5.24) 


Vn,l = 


i 


[1 - b{k)]^ {2Try 


(5.25) 


1 


( 2 d )2 


ln,l{0) — 2In,l{2ei) H- - — In,l{2e\ + 262 ) + ■ 


Bounds on the suprema of In,i{x), Kn,i{x), Tn,i{x), Un,i{x). In Section [331 we have 
bounded T-Ci’\x) in terms of SRW-integrals. To compute the bound on we need to rely on 

sup/„,z(a:), s\xi)Kn,i{x), sup Tn,iix), sup Un,i{x), (5.26) 

xgS xgS xsS xGS 

for different sets of vertices S. For hnite sets, we simply take the maximum of the elements. 
To obtain bounds for inhnite S, we use monotonicity of the SRW-integrals as formulated in 
the following lemma: 

Lemma 5.1 (Monotonicity of In,i{x) and L^^x) in x). Let n he a positive integer and 
consider x,y with xi > X 2 > ■ ■ ■ > Xd > 0 and 2/1 > 2/2 > • • • > 2 /d > 0. Then, 

In,i{x+ y) < In,i{x), Ln{x+ y) < Ln{x). (5.27) 
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This lemma is a combination of [271 Lemmas B.3 and B.4]. In the previous section, we 
have obtained bounds on Kn^i{x),Tn^i{x), Un^i{x) for a given x in terms of In,iix), Ln{x) and 
Vn^i- Thus, the supremum in (15.2611 is obtained for an x with the lowest order in the sense of 
Lemma l5.ll Therefore, when we use the following infinite set Q = {x € \xi\ > 2}, we 

bound the corresponding SRW-integrals by 

sup In,iix) = max{4^i(3ei), 4,i(2ei + 62 ), + 62 + 63 )}, (5.28) 

x€Q 

sup Lnix) = max{L„^i(3ei), Ln,i{2ei + 62 ), Ln,z(ei + 62 + 63 )}. (5.29) 

x&Q 


5.3 Bounds implied by the bootstrap functions and SRW integrals 

For the analysis presented in the previous sections, we use that bounds on the bootstrap 
functions /i, /2 and imply bounds on the NoBLE coefficients, see e.g.. Assumptions 12.71 
& 14.31 In the model-dependent papers, we prove that the coefficients can be bounded by 
combinations of simple diagrams. Simple diagrams are diagrams arising from combinations 
of two-point functions, like the triangle diagram that we have already seen for percolation: 


\ 7 p(®) — {GzGzGz)ix^• 


(5.30) 


The simple diagrams can then be bounded using the bootstrap function / 2 , e.g., 


Vp(^) 


' (—7r,7r)'^ 


Gzikfe^^-^ 


{2ttY 


< 




c{kY 


dYk 

{2ttY 


2d-2 
2d-I 



3 

73 , 0 ( 0 ), 


(5.31) 


which can be computed numerically. In this section, we present the bounds that we use in 
our implementation of the analysis. These bounds are optimized so as to obtain the best 
numerical bounds on (5 possible for Assumption 14.31 


5.3.1 Simple and repulsive diagrams 

In this section, we explain how to bound simple diagrams and repulsive diagrams, that we 
define below. 

The bounds on the simple diagrams are model dependent. However, for all models that 
we consider, the bounds are using the same idea that we present below. We require some 
notions that have not been introduced yet, namely, the minimal length of an interaction, the 
modified two-point function and repulsive diagrams. As an example, we give the definitions 
for percolation. The definitions for LT and LA are straightforward generalizations. 

For percolation, we define {x y} to be the event that there exists a path consisting of 
at least m open bonds connecting x and y. Similarly, we define {x y} as the event that 
there exists a path consisting of exactly m open bonds connecting x and y. [This is not the 
same as the event that the graph distance in the percolation cluster equals m, but it implies 
that it is at most m.] To characterize these interactions, we define the adapted two-point 
function, which are given for percolation by 

Gm,z{x) = (^x Y^ , Gm,z{x) = IPz (x y) ■ {5.32) 
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For all models under consideration the following holds, for all n G N, 

^ {2dflz)""cmix), (5.33) 

where Cm{x) is the number of m-step self-avoiding walks starting at the origin and ending at 
X. A diagram is repulsive if the paths involved do not intersect. For percolation, for example, 
the repulsive bubble and triangle are given by 

^rn„m2ix) = '^'^z ({0 v} O {y , 

^mi,m2,m3{x) = ^ (^{0 u} O {u O {y , 

where the symbol o denotes the disjoint occurrence, which is a standard notion in percolation 
theory, see e.g., [20]. For the examples above, it means that the occupied paths that are 
required to exist make use of disjoint sets of bonds. 

Below, we use the symbol (/ ® g){x) to describe that the paths involved for / are disjoint 
from the path involved in g. For example, {Gn,z ®Gm,z){x) represents the diagram where the 
path used in Gn,z is disjoint from that used in Gm,z- We define an{x) to be the number of n- 
step simple random walk from 0 to x that never use a bond twice, and note that {an®am){x) < 
an+m{x)- Further, {an®Gm,z){x) represents the combination of an n-step SRW path counted 
in an and a percolation path of length at most m, where, given the n-step SRW path, Gm,z{x) 
is the probability that an occupied percolation path exists that uses different bonds than the 
n-step SRW. The bounds that we explain below rely on the following bounds on the two-point 
function that hold for all models that we consider: 

Gm,z{,^) < ^m,0 < {2dfiz)'^{am ® Gz){x), (5.36) 

S—1 S—1 

Gm,zix) < '^Gi^zix) + Gs,zix) < '^{2di2zyai{x) + {IdflzYias IS' Gz){x), (5.37) 

i=m i=m 

^n,m{x) ^ {2dfiz) (®n ^ Gm,z){x) i^n+l,m{x), (5.38) 

for s, n, m G N with s > m. 


(5.34) 

(5.35) 


5.3.2 Bounds on simple diagrams 

Here we derive efficient numerical bounds on simple diagrams. Throughout this section, we 
fix n G N and ruj G N for i = 1,..., n. Further, we use the notation rriij = 

Assuming bounds on the bootstrap functions, we obtain the following bound for non- 
repulsive diagrams from (15.33^ : 


{Gm\ * Gm2 * ■ ■ ■ * Gr, 


< {2dllzr^’^{D*^^’-*G*^){x) 
2d 


< 


2d - 1 


2d - 1 


Ts K, 




(x). 


(5.39) 


When mi^n > 10, we also use this bound also for the repulsive diagrams. For mi^n < 10, 
we instead use the repulsiveness to reduce the numerical bounds in most cases by around 
50%. We obtain these improved bound by extracting short, explicit contributions. In high 
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dimensions, short connections typically give the leading contribution to diagrams, so treating 
them more precisely often pays off. This requires the computation of an{x), for all x € 
and n € N. For n < 10, we compute these values using a simple Java program that can be 
downloaded from the website of the first author |14| . 

We start by explaining the bound on the repulsive bound for the example of a bubble. 
We fix an M G N with M > mi + m 2 and use (15.381) to obtain 


M—m 2 —1 

—m 2 ,m 2 (^) 

i=m\ 

M—m 2 —1 M—l—si 

si=mi S2=m2 

+ {2dfl,)^{D*^*Gf){x) 

M-l 

= ^ {i + l-mi^2)aiix)fll 


(5.40) 


^=^1,2 


+ (M - mi,2){2dfl,)^{D*^ * G,){x) + {2dfl,y^^ * Gf){x). 


- \M / , /^'k2\ 


We extend this idea to obtain a bound on the triangle of the form 

M — 1 ^ — "^2,3 2 —m 3 —Si 

’^mi , m 2 ,m 3 ( 3 ^) ^ ^ ^ o,i{x)n^ E E 1 

2=mi,3 si=mi S2=m2 


(5.41) 


M—m2,3 —1 

+ ^ {M-m2,3-s){2dfl,)^{D*^*G,){x) 

s=mi 

+ (M - mi,3)(2dA^)^(T>*^ * Gf){x) + {2dfl^)^{D*^ * Gf )(x) 

(i + l-mi,3)(z + 2-mi,3) ^ 

= 2^ -^- ai[x)^i^ 

i=mi^S 

+ (M - mi^3){2dji,)^ * Gf)(x) + {2dji,)^*Gf){x). 


In the same way, we bound the square by 


mi,m2,m3,m4 ) (5.42) 

M-l 3 ..3 

- E ^II(^“”^h4 + s)«i(a;)M^ + -tl(M-mi,4 + s)(2(i/22)^(T)*^*G^)(x) 


6 ^ 

2=mi^4 s=l 

(M — mi^4)(M — 1 — mi^i) 


+ 


6 1 

S=1 

j- -.MfrM.M 


*Gf){x) 


+ (M - mi, 4 )( 2 dA^)^(Zl*^ *Gf )(x) + {2dJi,y^\D*^^ )(x). 


Jr, \M/r\*M 


5.3.3 Bounds on weighted diagrams 

Weighted diagrams, such as 1-C^’\x) in (I3.10p . are bounded using the bootstrap function f^, 
defined in (12.3p . It is especially beneficial to extract explicit contributions from the weighted 
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diagrams, as the bound produced by /s is not very sharp and decreases quite fast 

when we increase 1. We conclude from (I5.37P that, for I < M, 


M-l 

\\y\\lGMiG.^D*^){x -y)<Yl ( 5 . 43 ) 

y i=i 

For X that are close to the origin we can bound 'H^z''{x) quite efficiently. We abbreviate 
Hz{x) = ||x|| 2 G^(x) and compute that 

?7°’0(0) = 77,(0) = 0, (5.44) 

nY{Q) = {D^Hzm = ( 5 - 45 ) 

L 

77°’2(o) = (D*D* 77,)(0) = Y1 + e^), 

= — (4(2d — 2)G,(ei + 62 ) + 4G,(2ei)). (5.46) 

Most of the weighted diagram that arise in our bounds on the NoBLE coefficient are repulsive. 
Using this repulsiveness, we can make another substantial improvement. As an example, let 
us consider (oi 0 Hz){0). All connections counted in 77, need to make at least three steps as 
the direct step is used by oi, so that 

^(ai ® 77,)(0) = (77 ® 77,)(0) = G 3 ,.(ei), (5.47) 

which is numerically a factor l/2d better than the bound in (15.451) . If we consider two explicit 
steps of Ci + Ck, then we obtain the bound 

(aa ( 8 ) 77,)(0) <8d{{2d - 2)fll + G6,,(2ei)) (5.48) 

+ 8d{2d — 2)[fi^ + 4(2(7 — 3)/2^ + Cr 6 , 2 ( 6 i + 62 )). 

In our computations we have extracted all paths up to a length of six and manually computed 
the number of percolation paths that do not use any bond of the first path. This leads to 
excellent bounds on closed weighted repulsive diagrams. For bubbles and triangles, we extend 
the idea used in (j5.40p - (|5.4ip for mi, m 2 G N with mi + m 2 < 6 , to compute that 

5 

(Gmi,2 <8> Hz){0) < ^ plitti ® 772)(0) + (Ge,, ® 77,)(0), (5.49) 

i=mi 

5 

(Gmi ,2 < 8 ) Gm 2,2 ® 77,)(0) < ^ i(7 + 1 - mi - m 2)(7 + 2 - mi - m 2 )/i^(ai 0 77,)(0) 

i=mi-\-m2 

+ (6 - mi - m2)flliaQ 0 G, 0 772)(0) 

+/I®(a 6 0 G, 0 G, 0 772)(0). (5.50) 

The terms involving a* 0 77, are computed explicitly. All contributions involving a G, factor 
also contain a factor ae, numerically making them of order d~^. Since Hz{y) = \\y\\^Gz{y), we 
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can bound these minor contributions using /s by removing the added restrictions that the (8) 
convolution imposes and replacing it by a normal convolution. In this bound, we even bound 
o-eiy) < {2d)^D*^(y). Following this strategy, we reduce the effect of the bad bound /s and 
enhance our numerical bounds. 

5.4 Analysis of matrix-valued diagrammatic bounds 

In Section U we use the bounds in Assumption 14.31 to bound the terms appearing in the 
rewrite of the NoBLE equation. These bounds are stated in terms of the functions and 
2(N),q Bounds on these quantities are proved in the model-dependent articles, where they are 
stated in terms of matrix products, such as 

X 

iV-1 

^ < (N + 2)(h^B^w+ ^ , 

X M=0 

where h,v,w € and B,C € for some n >2 and A > 0. We need to sum these 

bounds over various sets of N to create and 

In this section, we explain how we compute the sum of these estimates. As a first step, 
we compute the eigensystem/spectrum of B, so the left eigenvectors fji and right eigenvectors 
(i to the eigenvalue A*. In our applications, there always exists a set of n independent left 
and right eigenvectors. As these vectors are linearly independent, there exists ri,... ,rn and 
bi,... ,bn such that 

n n 

v = '^riffi, w = '^biCi. (5.53) 

^=1 i=l 

We compute ri,..., using relations of the form 

n 

V = '^riffi = r]d, (5.54) 

i=l 

where = (ri,... ,rn), while the rth row of the matrix rj equals f]f. As the rows of the 
matrix y are independent vectors, rj is invertible, so that 

rj~^v = r, (5.55) 

which allows us to compute the Oj. The ftj’s are computed in the same way. We define 

Vi = Tiffi, Wi = biCi, (5.56) 

and note that these are also eigenvectors of B to the eigenvalue Aj. Thus, for A > 0, 

< vB^iv = (5.57) 

X i 

Using a geometric sum, we obtain 

^ E^^\x) < ^ =■ 

x,N X j * 


(5.51) 

(5.52) 
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To create a closed form expression for the bound on the weighted diagrams appearing in 
(15.521) and use that Vi and Wj are eigenvectors of B to obtain 

CXD CO 

N>0 X i N=0 i N=0 

CO A^—1 

+ E (^ + 2) E ( E ) C ( E ■ (5-59) 

N=0 M=0 i j 

For the second line, we rewrite the sums over N,M for fixed i,j as 

CO A^—1 CO CO 

£ (iV + 2) ^ CuJ, = £ £ (iV + M + 2)Af Af uf Cw,. (5.60) 

N=0 M=0 N=0M=0 

We use the geometric sum identity 

CO - 

E(’>+l)^” = 7rrAW. (5.61) 

n=0 ^ ' 

to bound (15.521) as 


El 

N,x 


|2^(n)/ 


(E E (^+2) + (E E (^ + 2) ) h 

i N={) i N=^ 

CO CO CO CO 

+ E E E (^ + l)Af Af FfCFJ, + ^ ^ ^ (M + l)Ar Af Ff CFJ, 

A^=0Af=0 ij N=0 M=0 


E'i’'A, 


+ E 


Wi 


1 


+ 


1 


(1-A,)2 ' 1-A, 
u?" Cw. 


:) +E aC^ 

1 


1 


+ 


1 


_ _L 

X n 


(1-A,)A1-A,) f-(l-\)(l-A,)2 

•'ij ^1j 


(1-Ai)2 ' 1-Ai 

_, Oabs 

Pas* 


(5.62) 


This highlight how we bound The bounds on are obtained in a similar way. 


Remark 5.2 (The numerics of the matrix powers and the eigensystem). In the classical 
lace expansion, similar exponential hounds appear for the Nth lace-expansion coefficient as 
a function of N, in which the base of the exponential is roughly bounded by the sum of the 
matrix elements J2ij (in fact, it is worse, since there loops do not have length at least 4). 
We use the matrix-valued bound to optimally use the fact that the number of steps in shared 
lines in loops appearing in our lace-expansion bounds provide information about the number 
of steps in other lines that are part of the same loop. Such bounds are most easily expressed 
in terms of matrix products. In these matrix bounds, the largest eigenvalue of B decides the 
magnitude of the bounds (see (I5.58p and (|5.62p ). For example, for percolation and d = 11, 
the matrix B equals 

/0.0134202 0.0112907 0.0257405\ 

B = 0.0127527 0.0108018 0.0338533 . 

\ 0.028009 0.0260537 0.0401418/ 
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For this matrix, the largest eigenvalue equals Ai ^ 0.073, which is quite small, certainly when 
compared to Yli j ~ 

6 Completion of the bootstrap argument and conclusions 

In this section, we complete the bootstrap argument, and explain how the conditions are 
verified in Mathematica notebooks. We start by summarizing where we stand. 

Verification of the bootstrap conditions. In Section [3l we have verified the bootstrap 
conditions. Part of this verification was the improvement of the bootstrap bounds on the 
functions /i, /2 and /a defined in (l2.1D - (j2.3p . the latter being technically the most demand¬ 
ing. A sufficient condition on the bounds on the NoBLE coefficients that allows us to improve 
the bounds on /i, /2 and /a is formulated in Definition 12.91 in terms of bounds on a simplified 
rewrite of the NoBLE formulated in (ll.37|) . that is formulated in Assumption 12.71 These 
bounds on the rewrite are reformulated in terms of the original NoBLE coefficients 
and Lf^’^ in Appendix [Dj The assumptions that we need to verify in the model-dependent pa¬ 
pers m and [18] are Assumptions 12.2112.41 and Assumptions 14.1114.31 The initial Assumptions 
I2.6H2.8I are replaced using Proposition 14.51 in Section HI 

Implications of the completed bootstrap and proof main result. The bound on 
/i implies a bound on Zc, while that on /2 implies that the infrared bound holds with an 
explicit estimate on the constant given by r2. Thus, /i and /2 imply our main results for the 
mo del-dependent analysis for percolation in [T8| and for LT and LA in m- The bound on /s 
implies bounds on simple weighted diagrams, such as weighted lines, bubbles and triangles. 
These bounds are crucial in improving the bound on / 2 , which implies our main result. 

We now discuss the improvement of /s in more detail, splitting between the model- 
independent and the model-dependent parts, and their relations to SRW integrals, the bounds 
on the NoBLE coefficients and the Mathematica notebooks that finally complete the analysis 
and thus complete our proofs. We begin by discussing the model-independent improvement 
of /a and SRW integrals. 

Model-independent improvement of /a and SRW integrals. In Section [3.31 we have 
proven that the bootstrap conditions on /a hold. Both for the verification of the conditions 
on the initial point in Section 13.3.31 as well as for the improvement of the bounds on /a in 
Sections [333H3331 we have formulated our conditions in terms of x-space SRW integrals such 
as In^i{x), Kn^i{x),Tn^i{x),Un^i{x), Ln{x),Vn{x) and Jn,i{x). The numerical values of these 
functions are crucial to allow us to verify that the required bounds on the initial point hold, 
and to improve the bound on /a. Thus, in order to perform a successful bootstrap analysis, 
we need to obtain rigorous numerical bounds on such SRW integrals. These numerical bounds 
are explained in detail in Section [S] using the ideas by Kara and Slade in [27] . These values 
are formulated in terms of integrals of Bessel-functions and are computed analytically, up to 
a specified precision in the Mathematica notebook available at [TTj. This notebook needs to 
be compiled before the model-dependent parts can be performed, as the model-dependent 
analyses use them. We will explain the Mathematica notebooks in the next paragraphs. We 
first explain how the bounds on /i, /2 and /s can be used to obtain sharp numerical bounds 
on the NoBLE coefficients. 
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Model-dependent improvement of /s and bonnds on NoBLE coefficients. Hav¬ 
ing accurate numerical bounds on all the SRW-integrals involved at hand, we can obtain all 
bounds on /s formulated in Section [3.31 see in particular (|3.87p . The remaining main ingredi¬ 
ent of the bounds on /s is formed by the bounds on the NoBLE coefficients. The initial bounds 
on fi, f 2 , /s imply bounds on Zc, Gz{k), as well as on several simple weighted diagrams implied 
by /s, in terms of ri,r 2 and Ts. These, in turn, allow us to prove bounds on the NoBLE 
coefficients and to identify for every > 0 and • G {H, S'", AH, {AH'", 0}, {AS'', 1}}, as 
formulated in Assumption 14.31 This is, next to the derivation of the NoBLE, the second main 
result in the model-dependent papers [18] and [T7], where [T 8 | treats percolation, and [T7] LT 
and LA. As a result, we then have all necessary bounds needed to verify the improvement 
of the bootstrap bounds. The numerical verification is performed in several Mathematica 
notebooks as we explain next. 

Model-dependent and compnter-assisted verifications nsing Mathematica. In or¬ 
der to complete the model-dependent and computer-assisted proof, we need to run the model- 
dependent Mathematica notebooks that can be found at the hrst author’s web page [Hj. We 
first need to choose the dimension in SRW_basic .nb and run the file. This creates the numer¬ 
ical input needed for the model-dependent NoBLE computations, which is done in the files 
Percolation.nb, LT.nb, and LA.nb. 0 When compiled, these model-dependent files compute 
whether P( 7 , T, •) holds for the given input, Ti, r 2 , Ta and c^, by going through the loop 
described in Figured] We first compute bounds on simple diagrams for the initial point zj, 
for which we do not need to rely on the bootstrap bounds in terms of ri,r 2 ,r 3 , but we can 
rely directly on the link to SRW as formulated in Assumption 12.21 For z G {zi,Zc), we do 
rely on the bootstrap bounds to conclude bounds on the bootstrap functions. When both the 
bounds for zj and for 2 G {zj, Zc), are bounded by the chosen Fj for i = 1, 2, 3, we conclude 
the existence of appropriate 7 *, e.g. 7 * = (Fj -|- (computed bound))/ 2 , we have proven that 
P( 7 , F, •) holds. This numerical verification completes the argument for the given model in 
the given dimension. In |18| and HZ], it is also explained how we can then use monotonicity 
in the dimension d to obtain the result for all dimensions larger than the specified dimension. 
Further, the model-dependent notebooks contains an algorithm that helps to choose the opti¬ 
mal value for the constants ri,r 2 ,r 3 and c^,Cn,/, 5 . This is explained in the implementation. 
See e.g., [TB] for a more extensive discussion for percolation. 

Complete version of the Mathematica SRW files. Next to the basic version of the 
SRW notebooks, we also provide a complete version, which should be used when dealing with 
the models in relatively low dimensions. The notebook SRW.nb serves two purposes. Firstly, 
it allows us to compute SRW integrals with the desired precision in dimensions close to the 
upper critical dimension. The file SRW_basic.nb uses built-in Mathematica functions, and 
can only be used for percolation in d > 15, as otherwise the desired numerical precision cannot 
be guaranteed. In the file SRW.nb, we use a Taylor approximation of the Bessel function, as 
explained in detail in m Appendix B], instead, that gives reliable results for d > 9. These 
computations make compiling SRW.nb take around an hour, while SRW_basic.nb is compiled 
in less than a minute. Second, SRW.nb allows the use of arbitrary index sets <S for f^, see 
(HAj). Using the basic version only the vertex sets 5 = 0 and S = 'L'^\ {0} can be considered. 

^At the time of publication the LT and LA Mathematica files are not made publicly available, as the 
corresponding article is not finalised. 
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These extensions are crucial to reduce the dimension above which our results apply. 
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Notation 

brief description 

defined in 

SRW 

simple random walk 

Section 1 1.2. 11 

NEW 

non-backtracking random walk 

Section 11.2.21 

D 

SRW step distribution 

CH) 

i, K 

direction of a bond l, k £ {±1, ■ ■ ■, ±d} 

Section |1.2.1| 

u,v,w,x,y 

points on the lattice: 


k 

Fourier argument, so k € (—7r,7r)'^ 

Section |1.2.1| 

z 

parameter of a generating function, z € [0, Zc] 


f*g,n 

convolution of functions 


dif,Af 

directional derivative and Laplace operator 

(13T2D 


SRW and NEW two-point functions 

ra. (fTTT^ 

J 

permutation matrix with entries (J)t,K = 

(11.151) 

D(A:) 

diagonal matrix with entries (D(/i:))^K = 

(11.151) 

X{z) 

susceptibility of a given model 

(USD 

G, 

two-point function of general model 

(fT:MD,(f02D, (1071) 

[I] 

[I] 

coefficient of the NoELE expansion 

(11.241), (liTl) 


coefficient of the NoELE expansion 

(11.251), g21) 


NoELE parameters directly connected to z 

(11.241),(11.251) 


numerator and denominator of Gz 

(11.331),(11.341) 

^^,Z 5 

SRW contributions in ^z{k) 

(11.351),(14.231), (14.241) 

(^F,Z1 Cf,z 

SRW contributions in Fz{k) 

(I1.36D, (I4.18D, (14.171) 

R<s>,z{k), Rp^zik) 

remainder of the split of ^zik) and Fz(k) 

(ID.6D. (11.361) 

fiif2, /a 

bootstrap functions 

(I2T])-(I23D 

111 

assumed/concluded bounds on the fi 

dUD-dOD 

/?. 

assumed bound on the coefficient 

Assumptions 12.714.31 

D^^^{k) 

EtiidsDik)? 

((301) 


SRW two-point function and integrals 

(I3.35I)-(I3.38I) 

-^(N) -r-r(0),t,K: 

mo del-dependent split of the NoELE coefficients 

Section 14.1.11 


Table 1: List of notation, that is used in at least two different sections. 
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B Proof of Lemmas l2.12[j27T3 

Proof of Lemma \2.12l Using a telescoping sum identity, we write 

1 - = 1 - + • • • + (B.l) 

We reorder the sum over x using the symmetry of x g{x), to obtain 

d /i—1 

^ g{x) [1 - cos(A: • x)] = ^ g{x) X] ( X] “ cos(A:^ • x^)] 

X X 11=1 j = l 

d 

X f_l=l 

We use 1 — cos(nt) < n^[l — cos(t)], see e.g. Lemma Ti. 131 below, and the symmetry of g to see 
that 

d 

^ 5 (x)[l-cos(A:-x)] < ^ 5 (x)^x^[l-cos(A:/,)] 


^l=l 




X X 


(B.3) 


□ 


Proof of Lemma \2.12[ We obtain (|2.25|) by taking the real part of the telescoping sum identity 

J i—1 

l-U* = ^[l-e'*»]]Je‘*U (B.4) 


2=1 


J = 1 


In the following, we use that | sin(x + y)\ < |sin(x)| + |sin(y)|, \ab\ < (a^ + 6^)/2 and 
1 — cos^(a) < 2[1 — cos(o)] to conclude from (|2.25ll that 


J 


J i—1 


1 — cos(t) < EH - cos(ti)] + EEi sin(tj)|| sin(fj)| 

2=1 2 = 2 j = l 

< ^[1 - cos(ti)] + ^ X] + sin2(tj)] 


i=l 

J 


i=2 j=l 
J 


i=l 


^[1 -cos(ti)] H-^^sin^(ti) < j^[l -cos(ti)]. (B.5) 

□ 


2=1 


2=1 


C Decomposition of AGz{k) for the bootstrap function /s 


We improve the bound on the bootstrap function /s in Section 13.3.51 using the decomposition 
of AGz{k) into the five pieces Hi{k) — H^{k) defined in p.52l) - p.56ll . as stated in (|3.57p . We 
create this decomposition by expanding AGz{k) in the form of the simplified rewrite, defined 
in ()1.35l) - ()1.37p . and group the terms with common factors. This decomposition is a key step 
in our analysis. The decomposition is long and tedious and, due to the many terms, prone to 
errors. For this reason, we explicitly perform it here in the appendix. 
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Expanding the terms. Our strategy is to first expand 

AG^{k) = +^z{k)( A-1- ) +2y" ds^,{k)ds( - 1 -V (C.l) 

l-F,ik) ^ l-F,{k)J ^ \i-F,ik)F 

and then group them in terms of their common factors. To not loose track of terms, we la¬ 
bel them using the line number in the formula in which they appear, e.g., (jC.bl al and (IC.Sl dl. 


Before stating this, we remark how we extract SRW-contributions C*{k) from inverse 
powers of (1 — Fz{k)) using the notation introduced in (|3.4Up - (l3.44p : 

^ ^ =C*{k)-E{k), (C.2) 


1- F^{k) 1 - F^(0) -k aF,^[l - D{k)] + Rf,z{ 0] k) 


1 


{l-Fz{k)r 

1 

(i-Fzik)r 


= C*{ky -E{k)C*{k) - 


E{k) 


= C*{kY - E{k)C*{kY - 


i-Pyky 

2 E{k)C*{k) 


E{k) 


i-Eyk) {i-Fz{k)y 


(C.3) 

(C.4) 


Further, we recall the form of ^z{k) and Fzik), see pi.35p - pi.36p . and note that 

A^z{k) = — a^^zD{k) + AR^^zik), AFz{k) =—ap^zD^k) + ARp^zik)- (C.5) 

Then, we proceed by expanding the terms in pC.ip . 


First term of (jC.ip . Using (IC.2p we compute 

A^zik) 


1-Ez{k) - 


= —ct<i,zF)ik)C*ik) + a^^zD{k)E{k) + 


—IcTla') 


:=ila6lb'i 


AR^^zik) 

1 - Fz{k) 
;=ESlc) 


(C.6) 


Second term of (IC.ip . We first note that 

1 _ AFz{k) 


A- 


i-Fz{k) {i-Fzik)y 


+ 2 


ELiidsFzik)y 

{i-Fz{k)y 


We use pC.3p . (IC.5h and ^z{k) = + a^^zF){k) + R^^z{k) for the term containing 

and obtain 


(1-F,(A:))2 ^'{1 - Fzik)y 

= -aF,zD{k){c^,z + a^,zD{k))C*{ky 

^ V 

:=l[cna'l 

+ aF,zD{k){c^,z + a^,zD{k)) ( E{k)C*{k) + ^ 

\ y-Fz[k))^ 


AFz{k) 

(C.7) 


: = lEl]b) 

Rt,zik)aF,zD{k) ARp^zik) 


+ 


{1-Fz{k)y (i-Fzik)) 


Gz{k). 


: = lEl]c) 


: = lEl]d) 
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For the term including the mixed derivatives we compute that 


2d,(fe) 


{l-MkW 


=2^z{k) 


Yfs=i{oiF,zdsD{k) + dsRF,z{k)f 


{l-F,{k)Y 


(C.8) 


‘2i(y,p “1“ “1“ 


EtidsPjkf 

(1-F,(fc))3 


+ 2^zik) 


Eti dsRF,z{k? + 2aF,,dsb{k)dsRF,z{k) 


{l-Fz{k)f 

= 2aF,,(c^,, + a^,zb{k))&^^{k)C*{kf 

: = lEla) 

+ a^^zF){k))b^-{k) IE{k)C*{kf + 




l-Fz{k) 


;=lEl]b) 

+ 2aF,,R^,zik)- - ^?-V7:::C*{k) -2«^,,<h,(A;)-- ^^j^E{k) 


{1-Fzik)y 


{l-Fz{k)Y 


: = El]c) 

+ 2aF,z{aF,z - 


:=lEl]d) 


b^'^^{k) 

a-Fz{k)r 


+ 2^zik) 


:=iUBe) 

Eti dsRF,z{k? + 2aF,zdsb{k)dsRF,z{k) 

il-Fzik)r 


■.=(UI\i) 

Third term of (IC.ip . The expand the last term of (|C.1D as follows 


d / 

s=l \ 


-Fz{k) 


O I o D + dsRF,z{k)) 

2 2^{oi<s>,zdsD{k) + dsR<s>,z\k)) - , 

^ [l-Ez(k)Y 


S=1 


= 2a^^zb^^'^{k)C*{kf -2a^^zb^'^’^{k)E{k)C*{k) -2a*,,T»"“(A:) 


m 


:=iUMa.) 

+ 2{aF,z - l)a^,zb^^Yk) 


:=lElb) 

1 

il-Uk)Y 


l-Ez{k) 


— ICilci 


;=(Eld) 

d 


+ 2 - 


^ ^ {dsRi.,z{k)aF,zdsb{k) + ds^z{k)dsRF,z{k)^ ■ (C.9) 


{l-Fz{k)Y^^ 


:=l[cnei 


Regrouping the terms. In this section we group the terms (IC.6[ al- (|C.9[ el to derive 

Hi{k)-H^{k). 
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Let us start by recalling that M*{k) = D{k) — 2D^^^{k)C*{k) and see that 


(ITIHl ai + dCHa) = -a^^,C*{k)M*{k), 

(IClTl ai + dCJla) = + a^,zD{k))C*{kfM*{k). 

Reviewing the definition of Hi{k) we see that 

(Qa) + dEIla) + dUSla) + (lEla) =-iLi(A:). (C.IO) 


Then, we proceed with the term H 2 {k): 

(iriHlbi + dUllb) = a^^,Eik)M*{k), 

mi^) + ^h) = a^,,{c^,, + a^,Mk))M*{k) IE{k)C*(k) +, 

V il-F,{k))J 

(ICTlci + raid = 

{l-F,{k)r 

raibi + dusib) + dUTib) + dusib) + (Eric) + dusic) = -H2{k). (c.ii) 

H^{k) captures the following terms: 

cad) + cae) , 

do + de) - ^<‘0 ■ 

raid) + (IUSle)+(IUdlc) + dUTld) = -H^{k). (C. 12 ) 


In H^^k), we collect the double derivatives of the remainder terms. 


rai d + raidi 


‘^F'S>,zik) /\Rp^z{k) ^ 

l-Fzik) {1-Fzik)) 


-mik). 


The two remaining terms consist of mixed derivatives and are put into H^{k): 

+ ^e) = -H^ik). 


(C.13) 


(C.14) 


This completes the derivation of the split into Hi{k) — H^{k). 


D Bound on the rewrite 

Here we prove Proposition l4.5H ii. Namely, we prove that the bounds stated in Assumption 
EH are implied by Assumptions I4.1TI4.31 by computing expressions for the bounds required in 
Assumption 12.71 one after the other. 

For better readability, we divide this into five steps: Step 1 contains the six relatively sim¬ 
ple bounds stated in (I2.8ll - (l2.10p . while Steps 2-5 show the more elaborate bounds stated in 
(I2.1ip - ()2.13D . The bounds involved are not difficult, but rather somewhat tedious. Through¬ 
out this section, we omit the subscripts 2 ; and write, e.g., fi = fiz and ft = flz- 
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Step 1: bounds stated in (I2.8p - (I2.10I) . 


a) The bound j3^ in (12.8p is also assumed in (|4.3UI) of Assumption 14.31 


b) We recall the definition of in (j4.23D . use the bounds of Assumption 14.31 and the fact 
that the coefficients are non-negative to conclude 


1 - /5Si, - < 1 + f>SSn + t=E^/941„- 


2d/i^ 


3(0) 




p C 

In this and the following lines we bound p by 2 d-i • 

c) For ap^z, defined in p4.18p . we use the bounds of Assumption 14.31 to obtain 


2d^ 


1 - 


(0) 


1-m2% 




2dfiz 

1 - /i 


1 + 7 + ^J‘|3• 


(1-0) 


1 


(0) 




We bound /r on the left-hand side of the inequality from below by /3 . 

d) The absolute value of a^^z defined in (14.241) is bounded by 

2 d //2 


z\ = max • 


+ T^<9gW.,.2<i/9£;b + (D.3) 




(D.l) 


(D.2) 


e) To bound Ylx k (x) we use that the coefficients are defined by an alternating series, as 
well as Assumption 14.2114.31 to obtain 


Y,nf\x) < <jiJ2 =(”’'‘( 3 ^) - 

X,K N=0 X,K X,N,K 

(DO 


N=0 


(D.4) 


f) In the same way, we obtain that 


E*;(^) 


N=0 



(D.5) 


Step 2: bounds stated in (12.IIP . In the steps 2-5 we create bounds for the remainder 
terms, as defined in Section 14.1.31 For the bounding it is useful to note the explicit formulae 
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for the remainder Rp^z{k) and R^^z{k)'- 


and 


-L W ,, „ 


n=2 ^ N=2 i 


+ 


1-/^2 


^(-l)A^(e-i^‘xhW’;,(A:) - fi¥ZAk)) 


Afelo,!} t 


fJ' 

(l-Ar 




'-0/1 


i-Oi'-l 


N=1 


(1 — /r2j2 \ 


1 - 0 , 1-1 




(1-^2)2 


(D.6) 


^^°{k) (e“‘*^^ofr°’^i(A;) - {k))) (e“'*^‘i - ^) 


^0,^1 


R.Ak)= E + + 

AfG{0,l} A^=2 n=l 

- T-^ E 'i'‘(t)(2‘(^e-“‘ - J-‘(t)) 

J- hz , 

oo 

^ N=1 

1 - /^z ^ 


(D.7) 


Now, we start by bounding and for n > 1. Starting from (14.7p . we take the absolute 
value into the sums and use Assumption 14.21 for all 4'^ and to obtain 


YAn{x)\<^lY Y. E 

X&A * l-0,";l-n Xi'. Ei^i=^ 


(4o,o + g|S(xo)|) 

\ — A 


1 — A 


n 


n—1 

X ( n(|2‘-{a:. + e..-.)l + (.|2-‘-‘(x^.)l)) 

S=1 

X (|S''"-i(x„+ e,„_i +e,„)| +^|H'-"-i(a:„ + e,„_J| 

+ ^J-A ''”'~^ixn + eA\ + AA 


(D.8) 


As we sum over all x, the sums factorize and we obtain 


E i^”Wi s 

irGZ^ 


1 + 


< 


1 - /i 

^(1+ iPt) 

1 fl 


Ei=wi)(t^) (n E I=‘-(..)|)E 

XT’ ■ n -1 xr* 1 .■ I 


S — 1 Xs As — l 


f 2dixl3p ^ ” 

V 1-^ 


(D.9) 
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In the same way, we obtain 


Using the geometric sum, we compute that 




V V IF (x)| < 1 + / 2d/2^g^ N2 

2^ 2^ I nl j| _ ^ _ 2dg/3-bs 1 _ j > 

n=2x&Zd ^1 -^ 

n=l,;gZd ^ V 


1-M 


1 -/r 


Using these bounds, we obtain that Rf,z (defined in (ID.bp i is bounded by 

/Qabs 

“ / za/// 

fj- 


Y- ,o . M / 2d/r l + pt f2dm^\^ 
^ |F.,,(x)|< i_^ iT^ 


cGZ"^ 


+ 


+ 


+ 


1 - 
2d^ 

2dfj, 

(TW 

{2dfj.)‘^p. 


i-M 


LXJ 

([ X] +^i^ + 13 b'') 


we{o,i} 


N=2 


%+2dfiY.Pi!^ 


N=1 


(2 + + 


{2d)‘^iJ? 

(1-Mf 


^abs ^abs 


and (given in (ID.7p i by 


E <f>S.+fiS. + E 4”’ + ^ 


N=2 


^ ^ 2dlipf- V 1-M 

J- 1-/. 


2dix 


2d^ 


+ (/3b4u + + (1 + 4’^^) • 


1 — ' 1 — 
This proves the bounds stated in p2.11h . 


N=1 


(D.IO) 


(D.ll) 

(D.12) 


(D.13) 


(D.14) 


Step 3: The first bound in (|2.12l) . We call a sum as appearing on the left hand side of 
P2.12I) a weighted sum, where the factor ||x ||2 is called the weight (see also Section [23]) . To 
obtain a bound on weighted sums, we split the weight using (12.2411 . We explain this now in 
detail for the following contribution to Yhx Il®ll 2 l-^*,. 2 (®)l- 

E 11“^ + y'& h'W(-(» - e.) - m2-‘(!/))| ■ (D.15) 

1 Uj 

x,y,L 

We use Assumption 14.21 to bound 'h(, by and obtain 

Ea < E 11^' + !/ll^|2(a^)l(|E‘(!/ -e.)\+ ^|E-‘(!/)l). (D.16) 

1 Uj 

' 'T* i; / 
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Then, we rewrite ||x ||2 using the equality in (|2.24l) as 

(IDJ^ < n Y1 11^1121^(3^)1 X](^ + /3)l^'(y)l 

i LX 


y,i- 


+ 




1 — 
2p, 


E l=WI E tellids'to - edi + m|2-‘(!/)I) 


y,i- 


+ -I _ 2 ''(^)l)' (^•^^) 

M • -1 

2=1 x 

Since H is totally rotationally symmetric, we know that X)x 3^*l^(3^)l ~ Thus, the last term 
cancels. We apply the bounds of Assumption 14.31 to obtain 

< -r^ ((1 + m)/31h/3S° + + M/31s‘,o)) • (D-18) 

i /r 

Now we extend this idea. First we use Assumption 14.21 to obtain 

/ _ X n+1 N 

-fJ- f 


^||x||i|$„(x)| <- 




Y1 II X]^*ll2(^o.^o + ^l^(3:^o)|) 


xo,---,Xn+i i=0 


n+1 


X n + (D.19) 

S = 1 ts-1 

Then, we split the weight using the equality in (12.241) . We note that x i-+ totally 

rotationally symmetric so that any contribution due to the second sum in the equality in 
()2.24p cancels. After the split of the weight, the sums factor and we obtain 


II l|2|^ / M ^ 

2^ ||x||2|$n(x)| < I 


+ (n + 1) 


^||x||i|“(x)|(^|HX 


n+1 


(D.20) 


t^)""‘e(^*«-+i=wi)(ei=‘(. 


['^\\x\\l{\'^''{x + e,)\+n\E '■{x)) 


< 


2d// 

1 - n 


L,X 

n+l 




Other contributions to (recall ()D.7D ) are bounded in a similar straightforward manner, 
leading to the final bound 


I|3^ll2l^<i>+(3^)l <^As,r + I^aIr + 


(D.21) 




n=l 

2djl 
I — 
2dn 


N=2 

n+1 


n + 1 


((1 + + +/?23+o)) 

2 _^2 {f^As--,R,i + ^ (/5L‘,. + ■ 


N=1 
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This proves (I2.12p with /Saa.* equal to the right-hand side of (ID.2ip . 


Step 4: The second bound in (|2.12p . We bound the weighted sum of in the same 
way as for R^^z- We require the following three additional bounds: 

E E + e.) <2d^ > 

E + e,) <{2dfli + 15^^'^) , 

X,L^K, 

E + e, + e,) <(24)2/2 ) . 

X,L,K 

Next we explain how to derive these bounds for the example of ()D.22h . 

First, we use (I4.29P and then (|2.24p to obtain 

E E + Ck) E Edl^ + + x^)E^^\x + e^). (D.25) 

xGl4^ ^ xGli^ ^ 

Then, we note that all terms containing a single cancel by the total rotational symmetry 
of H and S'", i.e., 

E -p e^) = E x^,E^^'>’‘'{x -p e^) = E + e, + e^) = 0. (D.26) 

X,K. X^K,,L X,hi,L 

Applying the bounds of Assumption 14.31 we obtain (|D.22I1 . 

For n > 2, we bound Fn as we have bounded $,1 in Step 3. Indeed, first, we use that 

n / _ \ n+1 

E ll^||ilc^(a^)l < E E IIE^*ll2(^ + i=(^o)i) (y^) 

x&'^ /-0,...,Ln X0,...,Xn 1 = 0 ^ \ /^ / 

n—1 

X (n(|2‘-‘0'. + e..-.)l+M|2-‘-‘(a^.))l) 

^ S = 1 

X +e,J\+fj.\E‘^"-^{xn + e,^_^)\ 

+ n\E'-’--^{xn + e,J\+fj.'^\E-‘-’--^{xn)\). (D.27) 

We split the weight using the equality in (I2.24p . use that the sums factor and then use the 
bounds stated above to obtain 

E ii^iiii'^xWi < (^y*' 

+ (n-l) (f + + /^/51=\o) 

/ 24/2 (Oa.b!>\ri-l 


(D.22) 

(D.23) 

(D.24) 
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The additional contribution Rp^z (recall (ID.61) 1 are bounded in a straightforward manner and 
we obtain the bound 






n=2 




(D.29) 


+ f^{n - 1) (f + + /^/31 b\o) 


n=2 

OO / „ , N n+1 /Oabs'vn-l 


-El 


2d^ 


+ 

+ 

+ 

+ 


n=2 

2d/i 

Ai 




(/3^ 


■[ E 

Afe{o,i} 


1 + /i 

(f + /3^“) + A^/SIbso + /3S“) 

CXD 

+ §E(^'!S + /5i''’ + M/s£’) 


l + M '■ 

(N) . o(N) 

A'P,b,/ ' A-'Ai'.h,// 


Af=2 


(l-/r2)2 

(2(i)^/i^/i 


OO 


Ar=l 


'.'^“7 ft /"^abs I flabs , «aba , , , flabs \ 

_ p2yi wAsSt ' Pas'-,o ' Ps'- ' PPab'-jO) 


«aba.ab. , i^df 


/Dabs /Dabs 

rPAsPH^ 


We complete Step 4 by defining /3 ah,f as the right-hand of ()D.29D 


-A /i/3Ai^o -A /3S=) . 


Step 5: The bounds in (|2.13|) . Now we compute a lower bound on Rp^zi^) — Rp^z{k)'- 

-Rf,2(0) — Rp^zik) = ^^Rp^zix)[l — cos{k ■ x)] (D.30) 

X 

which could also be negative. We could apply Lemma [2. 121 and bound (ID.dOp from below by 
—/3 ah,f[1 — D{k)]. Since this bound plays a central role, we prefer to use a better bounds. 
For this, we note that [1 — cos{k • x)] > 0 for all k and x. So we can create a lower bound by 
only consider the negative part of the function Rp^z- 

We are not able to identify the points x at which Rp^z{x) is negative. However, we can 
decompose all summands defining Rp^z in (ID.6p into a negative and a positive part as the 
coefficients were defined via alternating sequences. In this ways, we create the lower bound 
as follows: 


7T^EE(- 1)'^ rfl(N):''0,''l (0) _ g-K^lO+^'-llnt^A'-O/l 

“ A^ / 


N=11^0,Li 
> - 

> - 




(1-A72)2 

{2d)^^j2 


OO 

EE ]^(2N),to,tl _ g-i(Lo-|-fcii)n(2N),/0,0 


7V=1 tOi'-i 

2, - °° 


(1-Ar2)2 


E (Cv + 2d/3r')[l-^(fc)], 


N=1 
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(D.31) 






























where we applied Lemma 12.121 and (ID.24I) only in the very last step. We use this idea to 
bound all terms in (|D.6p . except the minor term Yl’n =2 ^n{k) and in this way obtain 


-Rf,z(0) — Rp^z{k) 

1 - D{k) 


< 


E 


n=2 




1 - /i 

+ ^(n- 1) 

2dfl 


(D.32) 


+ E 


+ 

+ 

+ 

+ 

+ 

+ 


n=2 

2dfj, 

I — Ijfl 


2d^ 

1 

(/3g^ 


n+l /oab3^r^-l 


(A 


-if + + /«/3AH^0) 


n=2 

OO / rs 1- \ n+1 //Qabs'i^—1 


if + Pt) + /3S^) 


(1-^2)2 

(2d/i)^/l 


1 — Ai/ l + /i ^ 

OO 

+ ^"3™.,.,,, + ? E (A”*" + + 3‘'3i”’) 

N=1 

/ OO \ 

®,„+{2dApx:(^s>.+4”’) 


Ar=l 


_ ^2)2 + /^ABSO + + l^/3ri“o) 

+ 1?) + 2/x/3rr;ds™“) 

/3°" + /r/3ri"o + + A^/^riv + /f^^/^Asso) 


(2d/i)^ 

(2d/A)^ 

(2d/A)^ 


/^even ^ ^even ^ + ^/ 3 °dd ^ ^ /U^/SIh^o) ' 


We modified the desired inequality (|2.13p here slightly for better readability. This proves the 
lower bound in (|2.13p with equal to the right-hand side of pD.32p . 
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